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ABSTRACT 

Blazars with strong emission lines were found to be associated mostly with broad-line type 1 Active Galactic 
Nuclei (AGN). Hitherto, evidence for blazars identified with Narrow Line Seyfert 1 (NLS1) AGN was limited to 
very few individual cases. Here we present a comprehensive study of a sample of 23 genuine radio-loud NLS1 
galaxies which have the radio-loudness parameters, the ratio of radio (21 cm) to optical (4400A) luminosity, 
greater than 100. The sample, drawn from the SDSS and FIRST, is homogeneous and the largest of this 
kind. A significant fraction of the sample objects show interesting radio to X-ray properties that are unusual to 
most of the previously known radio-loud NLS1 AGN, but are reminiscent of blazars. These include flat radio 
spectra, large amplitude flux and spectral variability, compact VLBI cores, very high brightness temperatures 
(10 11-14 K) derived from variability, enhanced optical emission in excess of the normal ionising continuum, flat 
X-ray spectra, and blazar-like SEDs. We interpret them as evidence for the postulated blazar nature of these 
very radio-loud NLS1 AGN, which might possess at least moderately relativistic jets. We suggest that those 
steep spectrum radio-loud NLS 1 AGN in the sample are of the same population but with their radio jets aligned 
at large angles to the lines-of-sight. Intrinsically, some of the objects have relatively low radio power and would 
have been classified as radio-intermediate AGN. 

The black hole masses, estimated from the broad Balmer line width and luminosity, are within 10 6 ~ 8 M Q , 
and the inferred Eddington ratios are around unity. Unless the black hole masses are largely under-estimated, 
our result stretches the low mass end of the black holes of luminous, fast accreting radio-loud AGN to a smaller 
mass regime (the order of 10 6 M Q ) in the black hole mass-radio-loudness space where other normal AGN 
are seldom found. The results imply that radio-loud AGN may be powered by black holes with moderate 
masses (~ 10 6 ~ 7 M Q ) accreting at high rates (Eddington ratios up to unity or higher). The host galaxies of a 
few nearby objects appear to be disk-like or merger; and some of the objects show imprints of young stellar 
populations in their SDSS spectra. We find that some of the objects, despite having strong emission lines, 
resemble high-energy peaked BL Lacs in their SED with the synchrotron component peaked at around the 
UV; such objects constitute an intriguingly high fraction of the sample. The radio sources of the sample are 
ubiquitously compact. They are smaller than at most several tens of kilo-parsecs, suggesting a possible link 
with compact steep-spectrum radio sources. Given the peculiarities of blazar-like NLS1 galaxies, questions 
arise as to whether they are plain downsizing extensions of normal radio-loud AGN, or whether they form a 
previously unrecognised population. 

Subject headings: galaxies: active — galaxies: Seyfert - galaxies: jets - quasars: general - X-rays: galaxies - 
radio continuum: galaxies 



1. INTRODUCTION 

As a minority of AGN, radio-loud (RL) AGN differ from 
their radio-quiet (RQ) counterparts mainly in possessing 
prominent radio jets and/or lobes that produce strong radio ra- 
diation. In observations, a widely used division between RL 
and RQ AGN is the radio-loudness parameter R «10, defined 
as the ratio between the radio 5 GHz to optical B-band lumi- 
nosity (Kellermann et al. 1989; Stocke et al. 1992), though 
it is still a controversy whether R has a bimodal distribution 
below and above this value. RL AGN are an important labora- 
tory to study the formation of relativistic jets, which is not yet 
understood so far (e.g. Blandford 2000; Celott & Blandford 
2001; Meier 2003). Nevertheless, it is widely thought that 
radio-loudness (formation of radio jets) is possibly depend- 
ing on the accretion rate/state (e.g. Ho 2002; Maccarone et 
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al. 2003; Greene et al. 2006; Kording et al. 2006), the spin of 
the black hole (e.g. Blandford 2000; Maraschi & Tavecchio 
2003; Sikora et al. 2007), black hole mass (M BH , e.g. Laor 
2000; Lacy 2001; Dunlop et al. 2003; Mclure & Jarvis 2004; 
Metcalf & Magliocchetti 2006, but see Woo & Urry 2002 for a 
different view), and host galaxy morphology (e.g. Balmaverde 
& Capetti 2006; Capetti & Balmaverde 2006; Sikora et al. 
2007), or a combination of some of these factors (e.g. Ballan- 
tyne 2007). 

As an important sub-class of RL AGN, blazar is a collec- 
tive term for BL Lac objects and flat-spectrum radio quasars 
(FSRQs). Blazars are characterised by flat radio spectra at 
above ^lGHz, fast variability, high and variable polariza- 
tion, superluminal motion, and high brightness temperatures 
(see e.g. Urry & Padovani 1995; Celott 2002, for reviews). 
They are now believed to be RL AGN with the orientation 
of relativistic jets close to the line-of-sight, and hence their 
non-thermal jet emission is highly Doppler boosted (Bland- 
ford & Rees 1978; Urry & Padovani 1995). Blazars have 
distinctive spectral energy distributions (SED) that are char- 
acterised by two broad humps in the log v — log vf v repre- 
sentation. The low-energy hump is commonly interpreted as 
synchrotron emission, while the high-energy one as inverse 
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Compton emission scattering off the same electron popula- 
tion which produces the synchrotron emission. Depending on 
the peak frequency of the synchrotron hump, BL Lac objects 
are divided into Low-energy-peaked BL Lacs (LBL) with the 
peak around IR-optical wavelengths and High-energy-peaked 
BL Lacs (HBL) with the peak around UV/soft X-ray energies 
(Padovani & Giommi 1995). Classical FSRQs have generally 
the peak frequencies of the synchrotron hump similar to or 
even lower than LBL. Interestingly, it was suggested that the 
whole blazar family can be arranged in a sequence from HBL 
to LBL and to FSRQ in the order of decreasing synchrotron 
peak frequencies and increasing source power — the so called 
blazar sequence 5 (Fossati et al. 1998). However, exceptions 
were also reported in recent years as the presence of powerful 
FSRQs with HBL-like SEDs (though this remains a matter of 
controversy), whose synchrotron hump is peaked around the 
UV/soft X-ray band (see e.g. Padovani et al. 2003, and refer- 
ences therein). Following Perlman et al. (1998), we term FS- 
RQs with HBL-like SED HFSRQs (High-energy-peaked FS- 
RQs), and those with LBL-like SED— the classical FSRQs— 
LFSRQs (Low-energy-peaked FSRQs), as parallels to HBLs 
and LBLs. 

Until recently, RL AGN (blazars) with strong broad emis- 
sion lines were found to ubiquitously have the Balmer line 
widths greater than 2000 km s" 1 in full width at half maximum 
(FWHM). It was found that the line width is correlated with 
the orientation of the radio axis, with narrower lines as well as 
stronger optical Fe II emission in objects with smaller viewing 
angles to the radio axis (e.g. Wills & Browne 1986; Jackson 
& Browne 1991). This was interpreted as arising from a disk- 
like emission line region. A marked absence of radio-loud 
quasars (RLQs) was noted below the 2000 km s" 1 line-width 
cutoff (Wills & Browne 1986), that is the characteristic line- 
width range for narrow line Seyfert 1 galaxies — a sub-class 
of type 1 AGN in contrast to the classical Broad Line AGN 
(BLAGN). The situation has changed since the detection of 
radio emission from NLS1 galaxies. 

Apart from narrow line-widths of the broad Balmer lines 
(originally defined as < 2000 km s" 1 , Osterbrock & Pogge 
1985; Goodrich 1989), NLS1 galaxies also show other ex- 
treme properties compared to normal broad line AGN (see 
Komossa 2008, for a recent review), such as strong permit- 
ted optical/UV Fe II emission lines (Boroson & Green 1992; 
Grupe et al. 1999; Veron-Cetty et al. 2001), steep soft X-ray 
spectra (Wang et al. 1996; Boiler et al. 1996; Grupe et al. 
1998), and rapid X-ray variability (Leighly 1999; Komossa 
& Meerschweinchen 2000). Observational evidence suggests 
that NLS 1 galaxies tend to have small black hole masses and 
high Eddington ratios (defined as the bolometric to the Ed- 
dington luminosity ratio), R E( \d = iboi/^Edd ~ 1 (e-g- Boro- 
son 2002; Collin & Kawaguchi 2004). In fact, they were 
found to locate at one extreme end, opposite to BLAGN, of 
eigenvector 1 of the correlation matrix which is believed to be 
driven primarily by R E dd (Boroson & Green 1992; Sulentic et 
al. 2000). NLS 1 galaxies were once thought to be radio-quiet. 
This was simply a consequence of the fact that NLS 1 galaxies 
have a low radio-loud fraction (Zhou et al. 2006; Komossa et 
al. 2006b) and the small number of NLS1 galaxies known till 
then. 

Previous studies of the radio properties of NLS 1 galaxies 

5 The blazar sequence can be explained as, for instance, the decrease of the 
energy of electrons emitting at the SED peaks with increasing energy density 
of the seed photons for Compton scattering (Ghisellini et al. 1998) 



are all based on small samples. These showed that the ra- 
dio sources, if detected, are compact (less than a few hun- 
dred parsec), at modest radio power (Ulvestad et al. 1995), 
of steep radio spectra (Moran et al. 2000), mostly in the RQ 
regime (Stepanian et al. 2003), and of low RL fraction (Zhou 
& Wang 2002). Until recently, there were only several RL 
NLS1 galaxies identified and studied individually (see Ko- 
mossa et al. 2006b, for a review and references therein). Stud- 
ies using (small) samples of RL NLS 1 galaxies were carried 
out by Komossa et al. (2006b) and Whalen et al. (2006) for 
non-radio selected and radio selected objects, respectively. 
These studies confirmed that, as in 'normal' NLS1 galaxies 
(their RQ counterparts), RL NLS1 AGN are accreting at a 
high rate close to the Eddington limit. More importantly, as 
pointed out by Komossa et al. (2006b), with relatively low 
M BH and high R, they occupy a previously sparsely populated 
region in the M^n-R diagram. Most of the sources of that 
sample are compact steep spectrum sources. It was confirmed 
that RL NLS1 galaxies are rare (~ 7%), and very radio-loud 
objects with R > 100 are even more sparse (see also Zhou & 
Wang 2002; Zhou et al. 2006), compared to about 10-15% for 
normal BLAGN and quasars (e.g. Ivezic et al. 2002). Several 
of these studies discussed the starburst versus AGN contri- 
bution to the radio emission, and concluded that it is AGN 
dominated. 

Of particular interest, a few outstanding RL NLS1 
galaxies came to be known to exhibit blazar properties. 
RXJ 16290+4007 (Schwope et al. 2000; Grupe et al. 2004) 
was known as a blazar and its X-ray emission is dominated 
by synchrotron emission (Padovani et al. 2002), or else by 
soft X-ray emission that is typical for NLS1 galaxies, per- 
haps from the accretion disk (Komossa et al. 2006b). The 
SED of PKS 2004-447 6 (Oshlack et al. 2001) was well mod- 
eled with the blazar type (Gallo et al. 2006). Using infor- 
mation on the radio flux variations, very high radio bright- 
ness temperatures (> 10 13 K) were inferred for J0948+0022 
(Zhou et al. 2003) and 0846+5 13 7 (Zhou et al. 2005)— both 
included in our sample, arguing for relativistic beaming (see 
also Wang et al. 2006). They also show some other blazar-like 
behavior. Recently, Doi et al. (2006, 2007) performed high- 
resolution VLBI observations for several radio-loud NLS1 
AGN (three are included in our sample 8 ) and found that they 
are unresolved with milli-arcsec resolutions, setting direct 
lower limits on the brightness temperatures in the range of 
10 7 - 10 9 K. They also found that inverted radio spectra are 
common in the radio-loudest objects. The authors suggested 
that Doppler beaming, presumably resulting from highly rel- 
ativistic jets, can explain naturally the observations. In this 
respect, the most remarkable object of this kind is perhaps 
2MASXJ0324+3410, found in our recent work (Zhou et al. 
2007). It showed rapid variability in the radio, optical, and X- 
ray bands, and was even claimed to be marginally detected in 
TeV 7-rays (Falcone et al. 2004). Its non-thermal SED clearly 
resembles that of HBLs. Moreover, 2MASX J0324+34 1 is 
hosted by a relatively small, apparent disk galaxy with one- 
armed spiral or ring galaxy morphology. These enigmatic 
objects present a challenge to current models of both NLS 1 
galaxies and blazars, and render a unique opportunity to study 
jet formation in black hole systems accreting at high rates. 

6 However, the NLS1 nature of PKS 2004-447 is not certain; see Zhou et 
al. (2003) for a short comment. 

7 A formerly known blazar; see Appendix B for details. 

8 J0948+0022, J1633+4718, and J1644+2619. 
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Apparently, they do not seem to belong to any types of AGN 
currently known. An immediate question is whether they are 
some exceptional individuals or there exists a population of 
such objects. What are their general properties in continuum 
and line emission, central engine and host galaxies? What are 
their relations with other types of AGN? 

Motivated by these issues, we compiled a RL NLS 1 galaxy 
sample from the Sloan Digital Sky Survey (SDSS) database, 
aimed at a systematic study with a large and homogeneous 
sample. Since, as the first step, we focus on genuine RL ob- 
jects, we consider in this paper objects with the radio-loudness 
parameter R > 100 (see § 2.2 for the definition of R) for the 
following reasons. Firstly, the classification of AGN into RL 
and RQ is somewhat ambiguous in the range R = 10-50 
due to the presence of (a small number of) so called radio- 
intermediate AGN spanning in between the bulks of the two 
classes, which may have a different origin from genuine RL 
objects (e.g. Falcke et al. 2006a). The same may also be true 
for NLS1 galaxies. Secondly, in R > 100 objects contamina- 
tion of radio emission from the host galaxy is negligible (see 
§6.2.1). Thirdly, their jet component, if present, can be de- 
tected and studied relatively easily in other wavebands than 
the radio, such as optical and X-rays. Last, but not the least, 
they are the least studied objects as a sample, due to their ex- 
treme rarity, in contrast to less RL NLS 1 galaxies as the bulk 
in the samples of Komossa et al. (2006b) and Whalen et al. 
(2006). We defer to a later paper a comprehensive and com- 
plete treatment of the SDSS-FIRST detected NLS1 galaxies 
sample, including radio-intermediate objects, to address ques- 
tions such as the RL-RQ dichotomy for NLS 1 galaxies. Our 
results show that there exists a population of NLS 1 galaxies 
that are blazar-like, despite their rarity. 

The compilation of the sample and the analyses of multi- 
waveband data are described in § 2. The broad band con- 
tinuum and emission line properties are presented in § 3 and 
§4, respectively. In §5 the properties of their host galax- 
ies are presented. We discuss other properties such as black 
hole mass and accretion rate, and isotropic radio emission, as 
well as the implications of the results in § 6. Throughout the 
paper, we assume a cosmology with Hq= 70 km s" 1 Mpc" 1 , 
Vt M = 0.3, and CI a = 0.7. Errors are quoted at the 68% confi- 
dence level unless mentioned otherwise. We use the following 
conventions to denote power-law spectral indices: a v in the 
frequency domain \S(v) oc v a " , specifically a r =ai„ in the ra- 
dio band], and a\ in the wavelength domain [S(X) oc X ax ]; 
and the X-ray photon index T, defined as f P ho(E) oc E~ r , 
where f p ho(E) is X-ray photon number flux density. 

2. SAMPLE COMPILATION AND MULTIWAVELENGTH DATA 
ANALYSIS 

2.1. SDSS NLS1 galaxy sample 

We have carried out a systematic search for NLS1 galax- 
ies from the SDSS spectroscopic samples and the first results 
have been published in Zhou et al. (2006). The procedures 
of data analysis have been documented in that paper in detail 
and only an outline is summarised here. We carefully fit emis- 
sion line spectra, AGN continua, and host galaxy starlight in 
a self-consistent manner. Firstly, galaxy starlight and AGN 
continuum, as well as the optical Fe II emission complex are 
modeled (see Appendix A for a brief account) and then sub- 
tracted. The optical Fe II multiplets are modeled with the Fe II 
spectral data given by Veron-Cetty et al. (2004) for both the 
broad and narrow components. Emission line spectra are fit- 
ted with the following models, using a code similar to that 



described in Dong et al. (2005). The Balmer emission lines 
are de-blended into a narrow and a broad component, which 
are modeled by a Gaussian and a Lorentz profile, respectively. 
All narrow emission lines are fitted with a single Gaussian ex- 
cept the [O III] AA4959/5007 doublet. Each line of the [O III] 
doublet is fitted with two Gaussians, one of which is used to 
account for a possible blue wing as seen in a few objects. As 
shown in Zhou et al. (2006), our analysis procedures yield 
reliable and accurate measurements of emission line parame- 
ters. 

The SDSS spectroscopic samples of both galaxies and 
QSOs in the redshift range of z < 0.8 were analyzed based 
on the SDSS data release 5 (DR5, Adelman-McCarthy et al. 
2007). Following Zhou et al. (2006), we classify NLS1 galax- 
ies as those having the "broad" component of H/3 or Ha which 
is detected at the 10cr or higher confidence level and is nar- 
rower than 2200 km s" 1 in FWHM. As discussed in Zhou et 
al. (2006), objects selected as such fulfill naturally the sec- 
ond of the conventional criteria of NLS1 classification, i.e. 
[O III]A5007/H/3 < 3, and the sample is reliable and uniform. 
As a result, about 3300 NLS1 galaxies are found, supersed- 
ing the previously published 201 1 drawn from the SDSS DR3 
(Zhou et al. 2006, hereafter Zhou06 sample). 

2.2. Search for radio counterparts 

Having compiled a large NLS1 galaxies sample from the 
SDSS, we search for their radio emission using data of the 
FIRST 9 survey (Becker et al. 1995). The survey has a spa- 
tial resolution of 5"and typical 90% source positional uncer- 
tainty ellipses less than 1". The identification is performed 
by matching the SDSS positions of NLS1 AGN against those 
of radio sources following the procedure used by Lu et al. 
(2007). This procedure was designed to search for SDSS 
quasars detected in the FIRST survey, and is summarised here. 
Since extended radio sources often have multiple (diffuse) 
components in morphology and sometimes their positions are 
vague to define, we treat extended sources differently from 
compact sources in search of radio counterparts. By compact 
FIRST source we mean that only one source is found within 
3'of the optical position and is unresolved; for them, we use 
a 2"matching radius. We search for possible extended radio 
sources associated with the NLS1 AGN in two steps. First, 
candidates are selected if they match one of the following cri- 
teria. (1) Only one resolved radio source is found within 3'of 
the optical position, and the optical position is located within 
the size of the radio source. (2) Two radio sources are lo- 
cated nearly symmetrically around the optical position within 
3'. (3) More than two radio sources are scattered around the 
optical position within 3'. Then we visually inspect 6' x 6' 
cutouts of the FIRST images centered at an object to reject 
false matches. As discussed in Lu et al. (2007), this approach 
has proved to be effective for finding extended radio sources 
associated with AGN; about 24% of the radio quasars in the 
sample of Lu et al. (2007) are resolved by the FIRST survey, 
among which half show FR II type morphology (Fanaroff & 
Riley 1974). 

This procedure yields a sample of SDSS NLS1 galax- 
ies detected in the FIRST survey. Interestingly, we note 
that all the radio sources associated with these NLS1 
AGN are compact in morphology. For each object we 
calculate the radio-loudness parameter defined as R\a= 
/ 1/ (1.4GHz)// 1/ (4400A), where the fluxes are in the rest 

9 Faint Images of the Radio Sky at Twenty-centimeters 
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frame of the objects 10 . The rest frame 4400 A flux densities 
are calculated from the SDSS g-magnitudes assuming a spec- 
tral slope av=-0.5. An index of -0.5 is also assumed for the 
K-correction for the radio fluxes. We focus on very radio-loud 
objects only with R\a> 100 in this paper, and will present the 
whole sample elsewhere. 

We note that three NLS1 galaxies with apparent Ria> 100 
stand out for their extremely red optical continua (a\> 2.7), 
namely, J095919.14+090659.4, Jl 1 1354.66+124439.0, and 
J233903.82-091221.3. Our original fit to the continuum spec- 
trum with a fixed power law slope a\ = —1.7 yields substan- 
tial reddening of E(B-V)=0.7-l. Among them, J233903.82- 
091221.3 (z=0.66) has been found to show convincing evi- 
dence for substantial dust extinction (Wang et al. 2005) 11 . 
Given the fact that the other two objects resemble closely 
J2339-0912 in their SDSS spectra, their optical light is most 
likely subject to heavy extinction as well. After correction 
for extinction, their radio-loudness parameters are largely re- 
duced to Ria— 10-20, below or merely close to the RL/RQ 
dividing line. We therefore exclude these three objects from 
our current sample. 

2.3. Very radio-loud NLS1 galaxy Sample 

The above selection procedures result in 23 NLS 1 galaxies 
having R\a> 100, which compose our sample. By exclud- 
ing less RL objects with R\a < 100 our sample can be re- 
garded as very radio-loud to some extent. We believe that the 
radio-loudness estimates of our sample are not significantly 
affected by optical extinction, as is discussed in § 6.2.2. The 
objects are listed in Table 1 , along with some relevant parame- 
ters of the continua and emission lines. As demonstration, we 
show in Figure 1 examples of the SDSS spectra and their best 
fit models for two representative objects; they are chosen to 
represent the ranges of some characterising parameters of the 
sample. Measurable starlight contribution is present in only 
one object, J1633+4718, accounting for at most ^20% of its 
observed spectral flux density at 5100 A. Most of our objects 
show strong Fe II emission complexes (see Section 4). The 
[O III] lines are weak in general, all having [O III]A5007/H/3 
< 1.2. All of the objects fulfill the conventional definition of 
NLS1 (Osterbrock & Pogge 1985; Goodrich 1989), and are 
therefore bona-fide NLS 1 galaxies. 

The redshift and radio-loudness distributions are plotted in 
Figure 2. Our sample has a relatively high redshift distribution 
peaked at z ~ 0.5 (the median). We calculate the Z?-band abso- 
lute magnitudes M B assuming ct\=—\.5 for the K-correction, 
which are listed in Table 1. The M B values range from -20.8 
to -25.6, with a median of -22.8. Among the sample, 10 
objects are as bright as M B < -23 and hence can be classi- 
fied as narrow-line type I quasars. Three objects in the sam- 
ple have been previously identified with RL NLS1 galaxies 
and studied in detail, namely, J094857.32+002225.5 (Zhou 
et al. 2003), J084957.98+5 10829.0 (Zhou et al. 2005), and 
J172206.03+565451.6 (Komossa et al. 2006a). In addition, 
several objects were known as general RL AGN and have 
been studied previously, whose properties are summarised in- 

10 It should be noted that our definition of radio-loudness, if 1.4, is related to 
the previously commonly used i?6 cm = /„(6cm)/f„(4400 A) (Kellermann et 
al. 1989) viaRi 4= 1.9R6cm- Therefore, our radio-loudness cutoff corresponds 
to «6cm > 50. 

1 1 In addition to the extremely red optical continuum, its SDSS spectrum 
shows signatures of reddening of broad emission lines, detection of strong 
Ca II and Mg II absorption lines of non-stellar origin (see Wang et al. 2005, 
for details). 



dividually in Appendix B. 

2.4. X-ray and UV data analysis 

We have searched for X-ray emission from the sample 
objects using ROSAT source catalogues of the RASS and 
pointed observations of both the RXP (Voges et al. 1999) and 
WGA (White et al. 1995) catalogues. The matching radii are 5 
times the given positional errors of ROSAT sources, which are 
typically 10"-20". For each match, we visually inspect their 
SDSS optical images to avoid spurious matches in which the 
X-ray source is actually associated with another nearby ob- 
ject within the matching radius. Of the 23 objects, 12 were 
detected in the RASS and 4 in pointed observations, among 
which 2 detected in both. Thus a total of 14 objects — more 
than half of the sample — were detected in X-rays with the 
ROSAT PSPC 12 . 

All of the objects except J1633+4718 do not have sufficient 
X-ray counts to allow proper spectral modeling. Following 
Schartel et al. (1996), Siebert (1996), and Yuan (1998), we 
estimate the X-ray photon index T from the two hardness ra- 
tios 13 , wherever available, assuming an absorbed power-law 
spectral model. Both T and absorption column density Nh 
can be estimated as free parameters, or the index alone as- 
suming Nh to be the Galactic column density. This method 
had been verified to be reliable and robust under the assump- 
tion that intrinsic spectral shapes are indeed a simple absorbed 
power-law (Brinkmann & Siebert 1994). The X-ray fluxes in 
the ROSAT PSPC band (0. 1-2.4 keV) are calculated from the 
count rates using the energy to counts conversion factor (ECF) 
for a power law spectrum and Galactic absorption (ROSAT 
AO-2 technical appendix, 1991). The ECF for each X-ray 
source is explicitly calculated from the ROSAT PSPC effec- 
tive area 14 , by using information on individual T if obtain- 
able, or using the mean otherwise (T = 2.49, see below). We 
list some of the information about the X-ray observations in 
Table 2. 

For objects not detected in X-rays, we determine upper lim- 
its on their X-ray fluxes from count rate limits set by the 
RASS. Since the X-ray backgrounds are very low in most of 
the RASS observations, detection of a source is mostly deter- 
mined by the source photon counts that follow Poisson statis- 
tics. We set an upper limit of source counts as 12 for an object 
not detected in the RASS, as often adopted (e.g. Yuan et al. 
1998). The corresponding count rate limit is then calculated 
by using the effective exposure time at the source position as 
available from the RASS exposure map. Then the flux limit 
for each object is calculated using the same method above and 
assuming T as the mean. 

Among the ROSAT detected objects, J 1644+26 19 was ob- 
served with Chandra with a large number of photon counts 
accumulated. Moreover, for J0849+5108 the hardness ratios 
(HRl = 0.97 ±0.11, HR2 = 0.28 ±0.15) yield an abnormally 
flat spectrum T = 0.6^ (see also Zhou et al. 2005). For 
this object we also investigate its ROSAT data in spite of the 
relatively low number of counts (~ 60). We perform X-ray 
spectral and timing analysis for J1633+4718, J1644+2619, 
and J0849+5108. We briefly describe here the data reduction 

12 Positional Sensitive Proportional Counter 

13 Defined as HRl = (C 2 -Ci)/(C 2 +Ci), ffi?2 = (C 4 -C 3 )/(C 4 +C 3 ), where 
C\, C 2 , C3 and C4 are the number of photon counts in the 0. 1-0.4 keV, 0.5- 
2.0 keV, 0.5-0.9 keV, and 0.9-2.0 keV bands, respectively. 

14 We used those that are appropriate for the RASS or pointed observa- 
tions, accordingly. 
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FIG. 1 . — Examples of the SDSS spectra of two radio-loud NLS1 galaxies in our sample. The left-hand side panels show the rest frame spectra and the fits 
to the continuum with a nuclear power-law (blue color) and optical Fe II emission (magenta) model, plus host galaxy starlight contribution (green) wherever 
non-negligible (see § 2.1 for a description of spectral analysis). The panels on the right-hand side show a close-up of the continuum-subtracted emission line 
spectra and the best fits in the H/3-[0 III] region. SDSS J1633+4718 (upper panel), a flat-spectrum radio source, has the lowest radio-loudness and redshift, and 
its spectrum is among the highest S/N; it has the largest fraction (~20%) of detectable host galaxy starlight contribution among the sample. SDSS J1047+4725 
(lower panel), a steep-spectrum radio source, has the highest radio-loudness and redshift, and its spectral S/N is among the lowest in the sample. 
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procedures and present results in § 3.3. We use version 12.2 
of XSPEC (Arnaud et al. 1996) for spectral modeling. 

J1633+4718 (RXJ 16333+4719): The target was observed 
with the ROSAT PSPC-b with an exposure time of 3732 sec- 
onds on July 24th, 1993 (Obs-ID: 701549; PI: N. Bade). We 
retrieved the ROSAT data from the archive and use the XS- 
ELECT (version 2.3) tool of FTOOLS (Blackburn 1995) for 
data analysis. The X-ray spectrum is extracted from a circle of 
150 "radius, and a background spectrum is extracted from an 
annulus with inner/outer radii of 195/300". There are <~ 976 
net source counts, and the count rate is 0.26 ±0.01 ctss" 1 . The 
spectrum in the 0. 1-2.4 keV band is binned to have at least 
30 counts in each bin. J1633+4718 was also detected in the 
RASS with a weighted exposure of 909 s, yielding ^ 185 net 
source counts and a count rate 0.20 ± 0.02 cts s" 1 . We extract 
the RASS spectrum following the procedure described in Bel- 
lonietal. (1994). 

J1644+2619: The object was targeted with the Chandra 
ACIS-S (PI: S. Laurent-Muehleisen) with an exposure time 
of 2946 s on June 02, 2003. The data were retrieved from 
the Chandra data archive and were reduced from the level-2 
data set following the standard procedure using CIAO (ver- 
sion 3.4). There are 553±44 net source counts detected, with 
a count rate of 0.188 ctss" 1 . The 0.2-5 keV spectrum (very 
few counts accumulated above 5keV) is binned to have at 
least 20 counts in each energy bin. 

J0849+5108: The object was observed with the ROSAT 
PSPC-b in a pointed observation as a target in 1993, April 
with an exposure time of 4496 s (Obs-ID: 700821). There are 
only 63 net source counts detected in the 0. 1-2.4 keV band. 
We fit the unbinned spectrum by minimizing the C-statistic 
(Cash 1979; Nousek & Shue 1989), which is valid for Poisson 
statistics appropriate to the low counts regime. 

We also searched for UV data from the GALEX (Martin 
et al. 2005) GR2/GR3 data release 15 . Twelves objects in our 
sample have available photometric measurements in the far 
and/or near UV imaging bands with an effective wavelength 
of 15 16 A and 2267 A, respectively. Their UV images are 
mostly point source like. The GALEX magnitudes (in the AB 
magnitude system) are corrected for Galactic extinction using 
A Fl]Y /E(B-V) = 8.376 and A Nl]y /E(B-V) = 8.741, following 
Wyder et al. (2005). 

3. BROAD BAND CONTINUUM RADIATION 

3.1. GHz radio emission 

3.1.1. Radio morphology and spectrum 

The FIRST radio images of our RL NLS 1 sources are ubiq- 
uitously unresolved at a resolution of 5.4", which is remark- 
able. This sets upper limits on the projected size of 10-40 kpc, 
depending on the redshift. Among our sample there are sev- 
eral objects observed with VLBI, including J0948+0022 (Doi 
et al. 2006), J1633+4718 and J1644+2619 (Doi et al. 2007), 
and J1505+0326 (Dallacasa et al. 1998), as reported in the lit- 
erature; they remain unresolved at resolutions ranging from 
a few milli-arcsec to a few tens of milli-arcsec, which cor- 
respond to several to several-tens parsec in physical scale. 
These observations set constraints on the brightness temper- 
atures of T B > 10 8 " 9 K (Doi et al. 2006, 2007). Since all the 
radio sources in our sample are compact, we use hereafter 
the average of the FIRST and NVSS 16 measurements as their 

15 http://galex.stsci.edu/GR2/ 

16 The NRAO VLA Sky Survey 
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1 .4 GHz flux densities in the following analysis; this is to min- 
imize the effects of fluctuations on the results caused by large 
flux variations. Similarly, at other frequencies, in cases where 
multiple epoch measurements are available at the same fre- 
quency, the average flux is used. 

Multi-frequency and multi-epoch radio data were collected 
from the NASA/IPAC Extragalactic Database (NED). Of the 
23 objects, 1 1 have flux measurements at 5 GHz, for which 
the 1.4GHz-5 GHz spectral indices a x are estimated, though 
the observations were not simultaneous. The 5 GHz fluxes 
and a r are listed in Table 1, and the histogram of a r is shown 
in Figure 3 17 . It can be seen that the radio spectra are system- 
atically flat, with a median of -0.30 and a mean of -0.13. 
Adopting the conventional dividing line between flat- and 
steep-spectrum sources, a x = -0.5, 8 out of the 1 1 objects with 
a T measurement turn out to be flat-spectrum radio sources 
(a r > -0.5). In particular, several objects show an inverted 
radio slope (a r > 0) near 5 GHz and toward higher frequen- 
cies based on simultaneous multi-frequency measurements, 
namely, J0948+0022 (Reich et al. 2000; Zhou et al. 2003; Doi 
et al. 2006), J1633+4718 (Neumann et al. 1994), J1644+2619 
(Doi et al. 2007), and J1505+0326 (Tinti et al. 2005). Only 3 
objects have a r steeper than, yet very close to, -0.5. 

Among the remaining 12 objects without available slopes 
near 5 GHz, 9 have lower frequency 327 MHz-1.4GHz in- 
dices measurable (see Section C), out of which 6 have flat 
spectra (o^m > — 0.5). Thus, it is likely that most of these 
12 objects may have the same flat (or even flatter) indices 
at higher frequencies as well, unless there is a break near 
1 GHz, as in Gigahertz-peaked spectrum sources (GPS). This 
fact supports further the flat spectral index distribution found 
above for our sample. Though in part of the sample the radio 
slopes are estimated from non-simultaneous data, the current 
a r distribution should be correct in the statistical sense; there- 
fore, we expect that the bulk of our sample objects should be 
flat-spectrum sources. As is known, radio sources with flat 
spectrum, which can be produced by relativistic jet models 
(e.g. Blandford & Konigl 1979), have a fairly good chance to 
be blazars (e.g. Fugmann 1988). 

Now we can elaborate the radio-loudness parameter cal- 
culation for our sample by using the averaged 1 .4 GHz flux 
densities and the estimated a r for individual sources for the 
K-correction (if not available, a r =-0.5 is assumed as above). 

17 For J0849+5108, the flux variations at 1.4 GHz were substantial and the 
spectral shape around 5 GHz was flat or even inverted (Arp et al. 1979); we 
thus use a r =0 for this object in the analysis below. 
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We also make use of the measured optical slopes as presented 
in § 3.2. The improved radio-loudness values are listed in Ta- 
ble 1 . The differences between the previous estimates and the 
improved values are small, typically within 10% in log/? 1.4; 
such differences are comparable with changes in R14 caused 
by source flux variations typical of our objects (see below). 

3.1.2. Radio variability and brightness temperature 

A number of objects showed significant flux variations 
based on multi-epoch observations taken at the same wave- 
lengths 18 . In Table 3 we list those having the significance of 
variations higher than ~ 3 a between two observations, where 

a is defined as a = {S\ - Si) / \J o^+of 2 . We also give in 

the table the amplitudes of variations A S, the fractional am- 
plitudes AS/ (S) where (S) is the average of the two fluxes, 
and the time spans Af between the two observational epochs. 
Their variability is consistent with the fact that the variables 
have ubiquitously flat or inverted radio spectra when avail- 
able. The fractional amplitudes given in Table 3 have a me- 
dian of 43%. In particular, several objects show large am- 
plitude variations (^40-75%) on relatively short timescales, 
J0948+0022 (reported by Zhou et al. 2003), J0849+5 108 (re- 
ported by Arp et al. 1979), and J 1505+0326 within a few 
years, and J1548+351 1 by 44% within 207 days (140 days in 
the source rest frame). 

Large amplitude variability on short timescales has been 
commonly used to set lower limits on the apparent bright- 
ness temperature (T B ) of a radio source, since the size of 
the variable part of the source can be constrained from the 
light-crossing time (e.g. Krolik 1999; Fanti et al. 1983; Gopal- 
Krishna et al. 1984). Assuming that the variable portion of the 
radio flux is emitted from a region with a size smaller than that 
corresponding to the light-crossing time, which is the time 
span between the two observations, the apparent brightness 
temperature is estimated to be 19 

> AP ve = 2D 2 AS, 
B ~ 27r 2 ^ 2 (Af) 2 {\+z)nkv 2 {At) 2 } 

where k is the Boltzmann constant, D L the luminosity dis- 
tance, AS„ the variable portion of the observed flux den- 
sity, AP ve the corresponding radio power at the emission fre- 
quency in the source rest frame, v the observing frequency 
and Af the time span in the observer's frame. The estimated 
Tb are listed in Table 3. 

Of the 8 objects, 6 have estimated exceeding the equipar- 
tition brightness temperature ~ 3 x 10 11 K as proposed by 
Readhead (1994), which is expected to conform by most 
radio sources. In particular, 4 of the objects, all of flat- 
spectrum, have Tb exceeding the inverse Compton limit ~ 
10 12 K (Kellermann & Pauliny-Toth 1969). The highest T B is 
found to be 10 14 K in J1505+0326, from two VLBI observa- 
tions at 8.3 GHz with the variability significance of 3a. The 

18 In cases where two observations have much different spatial resolutions, 
we consider the variability to be genuine only if the higher flux was measured 
at a higher resolution. This is to avoid detection of spurious variations caused 
by the possibility that a measured higher flux may be contributed by contami- 
nating extended/nearby emission if the observation is made at a lower spatial 
resolution. We take this approach even though this problem actually has little 
effect to our objects, all of which are compact on the scales concerned, e.g. 
several arc-seconds. Thus, our results are conservative regarding the number 
of variables, and variability amplitudes and timescales. 

19 We noted that there was an error in our previous calculation of Tb in 
Wang et al. (2006) and Zhou et al. (2003), which was also pointed out by 
Ghosh & Punsly (2007). We have corrected it in this paper. 



extremely high Tb values are commonly explained as emis- 
sion originating from relativistic jets (Blandford & Konigl 
1979; Jones & Burbidge 1973). For instance, taking the in- 
verse Compton limit as a conservative limit, the minimum 
Doppler factor can be estimated as 5^ = (T B / 10 12 K) 1 / 3 . This 
results in 5 m i n =1.5-4.7 for the 4 objects (Table 3), consistent 
with the range of the Doppler factors of relativistic jets in- 
ferred for classical RL AGN (Ghisellini et al. 1993; Jiang et 
al. 1998). 

3.1.3. Summary of the GHz radio properties 

(1) All the radio sources of our sample are compact, be- 
ing unresolved at the FIRST resolution of several arc-seconds. 
(2) The majority of objects with available radio indices are of 
flat-spectrum (around 1 .4-5 GHz), some even with inverted 
spectra. (3) Most of the flat-spectrum sources show signifi- 
cant variability, from which the brightness temperatures are 
inferred to be as high as T B > 10 11 K, and even exceeding the 
inverse Compton limit of 10 12 K in 4 objects. These high T B 
values may suggest the presence of at least mildly relativistic 
beaming of the radio emission. 

The properties of the radio emission at low frequencies of 
our sample objects are summarised in Appendix C, which 
is important for understanding the intrinsic radio power and 
radio-loudness of these objects, as discussed in Section 6.3. 

3.2. Optical continuum emission 

Our standard algorithm to select NLS1 galaxies as above 
does not give information about the optical continuum shape. 
However, for RL objects the optical continuum may devi- 
ate from the 'canonical' shape in consideration of a poten- 
tial contribution of jet emission. To properly take this effect 
into account, we parameterise the optical continuum shape es- 
pecially for our RL sample. We model the AGN continuum 
in the 2500A-5500A rest frame range with a power-law or a 
broken power-law, and repeat the above spectral fitting with 
the slope as a free parameter. The fits are performed in the 
same way as described in §2.1 assuming no intrinsic redden- 
ing. The improved fits give a more realistic description of 
the AGN continuum; while, on the other hand, we find that 
the new fits have almost no effects on the parameterisation of 
the emission line spectra. The optical continuum fluxes in the 
source rest frame are measured directly from the power-law 
fits of the AGN continuum. 

It has been established that the optical luminosity L opt of 
NLS 1 AGN is tightly correlated with the H/3 line luminosity 
Lh/3 based on large samples (Zhou et al. 2006); hence Lh/3 
is a good indicator of L opt . We test this scaling relation for 
RL NLS 1 galaxies as shown in Figure 4, where the measured 
5 100 A luminosities are plotted versus Lh/3- As a comparison 
the <~ 2000 predominantly RQ NLS1 galaxies in the Zhou' 06 
sample are over-plotted, along with the fitted AL\5ioo- L H/ 3 
relation (Zhou et al. 2006, their Eq. 5). As can be seen, the 
RL NLS 1 galaxies reveal the same trend of correlation, but, 
of particular interest, lie systematically above the known L opt - 
Lh/3 relation for normal NLS1 galaxies. We estimate the ratio 
of the observed luminosity to that expected from the H/3 lu- 
minosity (L A 510o/£a5ioo) usin 8 the A£-A5ioo- relation. We 
find that the ratios range from 0.93 to 5.2, with a mean of 
1.9. The ratios differ systematically between flat- and steep- 
spectrum sources as groups, with the mean being 2.4 for the 
former and 1.8 for the latter, and 1.6 for those without radio 
indices. 
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FIG. 4. — Relationship between the nuclear monochromatic luminosity at 
5100A and the H/3 luminosity for the radio-loud NLS1 AGN of our sample 
(filled dots: flat-spectrum radio sources; open circles: steep-spectrum radio 
sources; asterisks: no available radio indices). Over-plotted are NLS1 galax- 
ies in the Zhou' 06 sample (small dots), predominantly radio-quiet, as well as 
the best-fit relation (dashed). It shows that the radio-loud NLS1 galaxies have 
systematically higher observed luminosities in excess of that predicted from 
the H/3 luminosity of normal NLS1 galaxies. 

The optical continua in the rest frame 2500A-5500A range 
can be described with a power-law for most of the ob- 
jects. The fits are relatively poor for two objects, namely, 
J 1037+0036 and J 1634+4809, which show convex-shaped 
spectra in the optical-UV band with a drop short-ward of 
Mg IIA2800A in wavelength. For one object, Jl 138+3653, 
the spectrum cannot be fitted with a power-law due to a rapid 
drop in the blue part, which results possibly from mild extinc- 
tion and is discussed further in § 6.2.2. The best-fit slopes 
show a large range from modestly red (a\ =-0.2) to blue 
(a\ =-2.4) spectra. Their distribution is shown in Figure 5 
(excluding Jl 138+3653), which has a median of -1.54 and a 
mean of -1 .43. We compare this a\ distribution with that of a 
RQ NLS1 sample, which comprises 55 SDSS-selected objects 
with reliably measured optical slopes (dashed line in the fig- 
ure, Ai et al. 2008, for objects with multi-epoch data the av- 
eraged slopes are used). It should be noted that that RQ sam- 
ple has the redshift and luminosity distributions indistinguish- 
able from those of our sample [the two-sided Kolmogorov- 
Smirnov (K-S) test probability of 0.48 and 0.52, respectively]. 
It can be seen that, though there is a broad range of overlap, 
our very radio-loud NLS 1 AGN appear to have systematically 
bluer continua than the RQ objects. The RQ sample have a 
median of - 1 .24 and a mean of - 1 . 1 5 . The chance probability 
that the two distributions are the same is only P c hance^ 0.01, as 
given by the K-S test. A comparison with the smaller NLS 1 
galaxies sample of Constantin & Shields (2003), which have 
UV-optical slope measurements with the Hubble Space Tele- 
scope (HST), yields a similar conclusion; the closest value 
of their sample to ours is the continuum slope of the median 
composite spectra, a„~ -0.79 (ot\~ -1.21). 

3.3. X-ray radiation 
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FIG. 5. — Histogram of the power-law slopes of the optical continuum (rest 
frame 2500-5000A) for our radio-loud NLS1 sample (solid line; except for 
one object which is fitted with a broken power-law). As a comparison, the 
dashed line represents the histogram of ct\ for a radio-quiet NLS1 sample 
selected from the SDSS with redshift and luminosity distributions compatible 
to ours (from Ai et al. 2008). 

Assuming an absorbed power-law spectral model and using 
the aforementioned hardness ratio method, the effective pho- 
ton indices are estimated for most of the objects with valid 
ROSAT/PSPC hardness ratios extracted from the ROSAT 
source catalogues. The estimated absorption A^h values are 
found to be well consistent with the Galactic values within 
errors, indicating that there is no significant X-ray absorp- 
tion in these objects. Hence, we adopt the photon indices 
T (Table 2) obtained by assuming Nu =A^ al for better con- 
straining the parameters. Furthermore, for 5 of the objects no 
physically meaningful result can be obtained by this method, 
because either the spectrum is complex and highly deviating 
from the assumed absorbed power-law, or the measured hard- 
ness ratios are displaced outside the valid zone by Poisson 
fluctuations. For objects with multiple X-ray data, we use 
the weighted mean for T and fluxes in the following analy- 
ses. Results from detailed spectral analysis for the 3 objects , 
J1633+4718, J1644+2619, and J0849+5108, are summarised 
below. 

J1633+4718 (RXJ16333+4719): The spectrum deviates 
apparently from a simple absorbed power-law model, with 
prominent soft X-ray emission dominating E <0.5keV and 
a hard tail at E >0.5keV. The best fit is achieved using a 
model of a power-law plus a redshifted blackbody compo- 
nent, and the absorption Nu is found to be in excellent agree- 
ment with the Galactic A^ al =1-79 x 10 20 cm" 2 . Hence N H is 
fixed to A^ al . The data and the best fit model are shown in 
Figure 6. The rest frame blackbody temperature is found to 
be 32.5^g g eV and the photon index of the underlying power 
law r = 1 .37 ± 0.49 (90% confidence for one parameter of in- 
terest). The latter is less well constrained due to the limited 
bandpass of the PSPC in the 'hard' X-rays. No temporal vari- 
ability is found within the observation interval. The spectrum 
measured in the RASS has almost identical spectral shape as 
that of the pointed observation. The best-fit parameters are 
in excellent agreement with those obtained from the pointed 
observation, with T = 1 -47^ g2 and the blackbody tempera- 
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ture 3Ct[o eV (90% confidence level), though the errors are 
large given the small source counts. No variation in either the 
spectral shape or the flux is found between the RASS and the 
pointed observation. 

J1644+2619: A simple power-law model is found to be an 
excellent fit to the spectrum (see Figure 7), resulting in a re- 
duced x 2 =13 (for 23 d.o.f). The fitted absorption Nh is al- 
most identical to the Galactic value A^ al =5.12 x 10 20 cm" 2 . 
The best-fit photon index is steep, T = 2. 19 ± 0.27 (90% con- 
fidence for one parameter of interest). Adding an extra soft 
X-ray component (e.g. blackbody) does not improve the fit, 
though with two more free parameters; the power-law pho- 
ton index becomes flatter, T <~ 1 .8^3, but is subject to a large 
uncertainty range which encompasses the single power-law 
model. We therefore suggest that the simple power-law model 
with a relatively steep index is the most likely spectral shape, 
though models composed of a flatter power-law with a weak 
soft component cannot be ruled out from the current data. 
Compared to the observation in the RASS, the X-ray source 
was about 2.5 times fainter at the Chandra observation epoch, 
though the photon index seems not to vary. 

J0849+5108: The spectrum can be fitted with a power- 
law with absorption Nu close to the Galactic value A^ al 
=3.0 x 10 20 cm" 2 . Fixing column density at A^ al yields a pho- 
ton index T = 1 .77^' to (90% confidence for one parameter 
of interest). We note that the extremely hard hardness ratios 
given in the source catalogue are incorrect, which may re- 
sulted from inappropriate background subtraction around this 
object in the presence of nearby, seemingly extended X-ray 
emission. 

Figure 8 shows the distribution of the effective X-ray photon 
indices T, which span a wide range from 1.7 to 3.3. Since 
the uncertainties of most of the derived T are relatively large, 
we try to disentangle the intrinsic distribution of T from the 
measurement uncertainties, using the Maximum-Likelihood 
method introduced by Maccacaro et al. (1988). Assuming that 
the intrinsic distribution of V is Gaussian, we find that the 
distribution is intrinsically broad, with a standard deviation of 
0.39!°,;^ and a mean of 2.60!°,;^ (at the 90% confidence level 
for two interesting parameters). We conclude that the soft X- 
ray continuum shape, as described by the power-law photon 
index V, exhibits a considerable variety among our RL NLS 1 
galaxies, with both flat and steep spectra. 

It should be noted that the thus estimated T from hardness 
ratios are only an indicator of the overall spectral shape in 
the 0. 1-2.4 keV range. A steep T value does not rule out the 
presence of a flat 'hard' X-ray component if both a steep soft 
X-ray and a flat 'hard' X-ray component co-exist. As an ex- 
ample, for J1633+4718, a slope of T ~ 2.8 is derived from its 
hardness ratios if the spectrum is approximated as a power- 
law. This is particularly true for NLS 1 galaxies, since a soft 
X-ray excess component is commonly seen in them. Better 
quality X-ray spectra in a wider energy band than that of the 
ROSAT/PSPC are needed to determine the true underlying 
power-law in these objects. 

3.4. Broad band SED 

To quantify the broad band SED, we calculate the 
commonly used effective spectral indices between the 
5 GHz radio, 2500 A optical/UV, and 2keV X-ray band, 
namely, a IO , a ox , and a rx ; they are defined as a 12 = 
-log(5i,i/5 J/ 2)/log(^i/^2), where S v \ and S„2 are the fluxes 
at the frequencies v\ and ^2, respectively, in the object rest 



frame 20 . The rest frame 2500 A fluxes are computed from the 
M-band PSF magnitudes, which have the effective wavelength 
(3542 A) very close to the redshifted 2500 A light (within a 
range of 500 A); for the K-correction, the measured optical 
slopes for individual objects are adopted (for one object with- 
out the slope available, ot\ =-1.5 is assumed). For the 3 ob- 
jects for which X-ray spectral fits were enabled, the X-ray flux 
densities at 2 ke V (/acev) are measured directly from the best- 
fit spectral model; while for the others /2kev are computed 
from the 0. 1-2.4 keV integrated fluxes. In doing so, as well 
as the K-correction, the derived T values are used for individ- 
ual sources whenever available, or the mean (T) = 2.6 is used 
otherwise. For those without X-ray detection, lower limits on 
a ox and a IX are calculated using the upper limits on the X-ray 
fluxes constrained by the RASS. The obtained a ro , a IX , and 
a ox are listed in Table 4. 

In Figure 9 we plot our RL NLS 1 galaxies on the commonly 
used a ro -a ox blazar diagnostic diagram, separating them into 
three groups: flat (underlying) spectra (V < 2.0, filled dots), 
steep spectra (r > 2.0, open circles), and those without V 
estimates (asterisks). For comparisons, we over-plot RLQs 
from the ROSAT detected quasar sample of Brinkmann et al. 
(1997). That sample is composed of both FSRQ and steep- 
spectrum radio quasars (SSRQ), whose loci largely overlap 
with each other. It can be seen that some of the objects do fol- 
low the trend defined by FSRQ/SSRQs, indicating their SEDs 
are similar to those of FSRQ/SSRQs. Specifically, all the 3 
objects in our sample known to have a steep radio spectrum 
are well consistent with normal SSRQs. It is known that FS- 
RQs (and SSRQs) occupy almost the same region on the a ro - 
a ox diagram as LBLs, owing to their similar SEDs (Padovani 
et al. 1997a; Brinkmann et al. 1996; Sambruna et al. 1996; 
Fossati et al. 1998). 

However, a considerable fraction of the objects clearly off- 
sets. We expect that these outliers must have different SEDs 
from those of classical FSRQ and LBLs. By invoking the 
classification scheme of BL Lac objects, we divide the a I0 - 
a ox diagram into the loci of LBLs and HBLs. A simple yet 
commonly used criterion separating HBLs and LBLs is a rx ; 
a nominal value is set to be around a' rx <~ 0.75 - 0.78 (de- 
fined between 5 GHz and 1 keV) suggested by various au- 
thors, with HBLs having smaller (flatter) a rx values. We adopt 
0.78 as the dividing line (as in, e.g. Padovani et al. 2003) , 
which corresponds to a rx =0.787 for the a rx definition adopted 
in this paper and is represented by the dashed line in Fig- 
ure 9. Following Padovani et al. (2003), a more elaborated 
locus for HBL is also indicated schematically, which is pop- 
ulated by X-ray selected BL Lac objects (e.g. Brinkmann et 
al. 1996). Of particular interest, most of the objects offsetting 
from the classical FSRQ and LBL locus actually fall within 
the HBL locus 21 . This indicates that the broad band SEDs 
of these NLS1 AGN are similar to those of HBLs. Follow- 

20 There are slight differences in the frequencies adopted in the a I0 and 
a ox definitions among the literature. Whenever we cite a IO and a ox which 
have different definitions, we first convert their values to those according to 
the definition used in this paper. 

21 In a strict sense, the HBL locus in Figure 9 represent the SEDs of the 
'pure blazar' component only or dominated by the non-thermal continuum. 
Our NLS 1 objects are, however, strong-lined and a significant contribution 
is expected from the thermal ionizing continuum to the observed optical flux 
(~ 50% on average, see § 3.2). For rigorous comparisons with HBLs, we 
should consider the 'corrected SED' corresponding to the pure non-thermal 
continuum only. However, such corrections are very small (Aq ox ~ -0.12 
and Aa ro ~ 0.05 on average) for the bulk of the objects, and the results are 
essentially not affected. 
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FIG. 6. — X-ray spectral fit for J1633+4718. Left panel: ROSAT PSPC spectrum, the folded best-fit model and the residuals; right panel: best fit model of a 
power-law plus a black-body (with Galactic absorption). 
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FIG. 7.— Chandra ACIS X-ray spectrum of J1644+2619 and the best- 
fit model as a steep power-law (r = 2.19 ± 0.27, 90% errors) with Galactic 
absorption. 



ing our terminology introduced in § 1, we refer to them as 
NLS1 HFSRQs, and those with SEDs typical of classical FS- 
RQs (LBL-like) as NLS1 LFSRQs. These NLS1 HFSRQs 
have ubiquitously flat radio spectra whenever available, con- 
sistent with a potential blazar nature. Those without X-ray de- 
tection (arrows) have their a IO and a ox limits consistent with 
the a vo and a ox values for either LFSRQs/SSRQs or HFS- 
RQs. Also, a few objects lie very close to the border line, 
on both the HBL and the LBL sides, suggesting a smooth 
transition between the two types. As a comparison, we also 
mark in the plot the position of 2MASX J0324+34 1 0, an ex- 
treme NLSl-blazar composite object, which has a HBL-type 
SED with the synchrotron peak in the UV/soft X-rays (Zhou 
et al. 2007). Clearly, 2MASX J0324+34 1 fulfills the HBL 
classification in terms of its radio-optical-X-ray SED, as ex- 
pected. We note there are a few objects lying in the vicinity 
of 2MASX J0324+3410 in the a ro -a ox diagram, suggesting 
their similarities and thus the previously known NLSl-blazar 
hybrid 2MASX J0324+34 1 is not unique. 

FSRQs/LBLs and HBLs are known to differ in their X-ray 
spectral shapes as FSRQs/LBLs have relatively flat spectra 
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FIG. 8. — Distribution of the effective photon indices of representing power- 
law of the X-ray spectra of our radio-loud NLS 1 galaxies, where available. 
All the indices are measured with the ROSAT PSPC in 0. 1-2.4 keV, except for 
J1644+2619 with the Chandra ACIS in 0.2-5 keV. For J1633+4718 we adopt 
T = 2.8 obtained from the hardness ratios to represent the overall spectral 
shape. 

while HBLs have steep spectra (e.g. Worrall & Wilkes 1990; 
Brinkmann et al. 1996; Padovani & Giommi 1995; Wolter et 
al. 1998); specifically, Y ~ 2.0 were found for FSRQs/LBLs 
and r ~ 2.5 or steeper for HBLs in the ROSAT band. As such, 
the overall spectral shapes from the optical to the X-ray band 
are concave in FSRQs/LBLs (a ox > a x , where a x = T-l is 
the energy index), while they are convex (a ox < a x ) in HBLs 
(e.g. Sambruna et al. 1996; Padovani et al. 1997a; Perlman et 
al. 2005). Parallelly, the X-ray emission of HFSRQs is hy- 
pothesized to be dominated by synchrotron radiation with a 
steep X-ray spectrum, as claimed to be seen in a few individ- 
uals (Padovani et al. 2002). For the NLS1 HFSRQs in our 
sample with measured X-ray slopes, we find that all but one 
(J1633+4718) have steep X-ray spectra (T > 2.0). We also 
examine the optical-to-X-ray spectral shapes for the objects 
with available a x by plotting a x versus a ox in Figure 10. As 
can be seen, the two objects (J1633+4718 and J0849+5108) 
with flat X-ray spectra a x < 1.0 show prominently concave 
optical-to-X-ray spectral shapes; whereas most of the steep 
X-ray spectrum objects are consistent with convex spectra, in- 
cluding 2MASX J0324+3410 with T ~ 2.2 (Zhou et al. 2007). 
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FIG. 9. — a ro versus a ox for the radio-loud NLS1 sample; filled dots: 
objects with flat X-ray spectra (r < 2); open circles: objects with steep X- 
ray spectra (T > 2); asterisks: objects without available X-ray slopes. Ar- 
rows represent lower limits of a x for X-ray non-detections. For compari- 
son, we over-plot ROSAT detected radio-loud quasars from Brinkmann et al. 
(1997), consisting of both FSRQs (small crosses) and SSRQs (small dots). 
The blazar-NLSl hybrid object 2MASX J0324+34 1 discovered by Zhou et 
al. (2007) is also indicated for comparison (square). The dashed line repre- 
sents a constant a IX =0.78, which is commonly used to divide BL Lac objects 
into LBLs and HBLs (see text for details). A more conservative, schematic 
locus for HBL is also indicated following Padovani et al. (2003, solid line). 
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FIG. 10. — X-ray spectral energy index ax (T~ 1) versus a ox . The overall 
optical-to-X-ray spectral shapes are divided into convex shapes (a ox < ax) 
and concave shapes (a ox > ax) by the dashed line (a ox = ax)- The numbers 
mark the object IDs as given in Table 1. 2MASX J0324+3410 is also plotted 
(square). 

We conclude that some of our radio-loud NLS1 AGN have 
SEDs similar to those of LBLs/RLQs, while the remaining 
similar to HBLs/HFSRQs. The X-ray spectral shapes are also 
broadly consistent with these two types of SED. 

Finally, in Figure 1 1 we plot the broad band SED in the 
\ogv - \ogvL v representation for the 14 objects with X-ray 
detection. Though there is only one broad radiation hump 



detected (due to a lack of hard X-ray and 7-ray data), they 
do reveal some similarities with the low-energy (synchrotron) 
hump characteristic of blazars. In particular, some of the 
HBL-like candidates are found to have the radiation hump 
peaked around the UV, supporting the results from the above 
a m -a ox distribution; these include Jl 146+3236, J1238+3942, 
J 1644+26 19, and J 1722+5654, etc. They also show SEDs re- 
sembling that of 2MASX J0324+34 10. 

4. EMISSION LINES AND AGN EIGENVECTORS 

The Eigenvector- 1 space involves primarily a set of well 
known correlations among the following parameters: the line 
width FWHM(UP), R 451 o, [O III]A5007/H/5, and soft X-ray 
spectral index (Boroson & Green 1992). ^4570 is the relative 
strength of the Fe II multiplets in terms of the Fe II-to-H/3 flux 
ratio, /?4S7o = Fe II/H/3, where Fe II is the flux of the Fe II mul- 
tiples in the range of 4434^-684A, and H/3 is the total H/3 flux 
(Veron-Cetty et al. 2001). [O IIIJA5007/H/3 is the strength of 
[O III] relative to that of H/3. Here we examine these cor- 
relations for our RL NLS 1 galaxies, and compare them with 
those for the large NLS1 sample of Zhou'06, which are pre- 
dominantly radio-quiet (plotted as small dots in the figures in 
this section). 

Figure 12 shows the correlation between the H/3 line- width 
and the relative strength of Fe II; again, a significant anti- 
correlation is present for the RQ NLS 1 sample (with a large 
scatter though), as shown in Zhou et al. (2006), whereas no 
significant correlation is found for the RL objects. A sim- 
ilar result is found for the correlation between ^4570 and 
[O ffl]A5007/H/3, as shown in Figure 13. While an anti- 
correlation between the Fe II and [O III] strength is evident 
for 'normal' NLS1 galaxies, there appears no such correlation 
for the RL objects. However, it should be noted that, given the 
weakness of most of these correlations found for the RQ sam- 
ple, our RL NLS1 sample is too small to draw statistically 
significant inferences. Also, we note that the high ^4570 val- 
ues for the RL objects span a relatively narrow range of ^4570 
in the diagram, which may partly account for the lack of the 
#4570-[O III]A5007/H/3 correlation. 

Another well known Eigenvector- 1 correlation is the anti- 
correlation between soft X-ray photon index T and H/3 line- 
width (e.g. Wang et al. 1996; Boiler et al. 1996), that was 
found to extend down to FWHM ~ 1000 kms" 1 (Zhou et al. 
2006, their Figure 17). We examine this property for the 14 
RL NLS1 galaxies with measured T, and find that they do 
not follow the T-FWHM trend defined by their RQ counter- 
parts. However, the sample size is small, and the scatter in 
the known relation is huge. Alternatively, the finding could 
be attributed to a different X-ray radiation mechanism domi- 
nating the ROSAT PSPC band in RL and RQ NLS1 galaxies 
(see below). 

We also search for potential correlations between radio- 
loudness R\a and ^4570 as well as [O IIi]A5007/H/3 within 
our RL NLS 1 sample; no significant correlations are found, 
however. 

In the process of optical spectral analysis, we have no- 
ticed that the Fe II multiplets are generally strong in our 
RL NLS1 galaxies. We plot in Figure 14 the distribution 
of /?457o, and, as a comparison, also the ^4570 distribution 
for the NLS1 galaxies in the Zhou'06 sample (normalized 
to the number of RL objects). The two-sided K-S test test- 
ing if the two distributions are drawn from the same popula- 
tion yields a probability level of 0.2% only; it is evident that 
the RL NLS 1 galaxies have systematically higher ^4570 val- 
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FIG. 1 1 . — Broad-band spectral energy distribution in the log vL v — logf representation for radio-loud NLS1 AGN with available X-ray data. The luminosities 
are at the emitting frequencies in the rest frame of the objects. The radio and infrared data are collected from the NED and supplementarily from recent radio 
surveys such as FIRST. The optical data are SDSS PSF magnitudes in the five bands (u,g,r,i,z). The UV data are GALEX near- and far-UV magnitudes. The 
optical and UV magnitudes have been corrected for Galactic extinction. The host galaxy contamination should be negligible, as shown in Section 2.3. For objects 
whose X-ray spectrum can be represented by a power-law, the estimated slope and uncertainties are indicated as a bow-tie. For J1633+4718, the X-ray spectrum 
is complex and the unfolded spectrum is presented. 

ues than their RQ counterparts. The average ^4570 is 1.03, 
compared to 0.82 for the 'normal' NLS1 galaxies with reli- 
able Fe II measurements in the Zhou' 06 sample. To exam- 
ine whether the large R4570 values are caused by systemat- 
ically weak (total) H/3 emission, we compare the H/3 flux 
distribution of the RL and RQ samples; we find that the 
two H/3 flux distributions are indistinguishable between the 
two samples (the K-S test). Moreover, the RL NLS1 sample 
has FWHM(Hf3) and [O IH]A5007/H/3 distributions indistin- 
guishable from those of the Zhou'06 sample (the K-S test). 
This result, together with the absence of the anti-correlations 
of /?457o-FWHM(H/3) and fl 4 5 7 o-[0 III]A5007/H/3 indicate 
that the high ^4570 values for the RL sample are real, rather 
than spuriously induced from dependence on other parame- 
ters. We thus conclude that the optical Fe II emission in our 
very RL NLS 1 galaxies is on average stronger than that in RQ 
NLS1 galaxies. In fact, this effect has been noted in the sam- 



ple of Komossa et al. (2006b, their Figure 3), but the statistical 
significance is not as high due to a factor of 2 smaller of the 
sample size. 

We have identified four [O III] blue outliers in our sample, 
with the [O III] lines significantly blueshifted by V([0 III])> 
200 km s" 1 relative to [Q II] 22 . The phenomenon of blue out- 
liers has been observed in a few radio galaxies (e.g. Tadhunter 
et al. 2001), and is especially frequent among NLS1 galax- 
ies (Komossa et al. 2007, and references therein). We find a 
fraction of blue outliers with V([0 III])> 200 kms" 1 of 17% 
among our sample, about a factor of 2.4 higher than that in 
the RQ NLS1 sample of Komossa et al. (2007). These au- 
thors also discovered that, for RQ blue outliers, the velocity 
offsets are correlated with ionization potential and line width. 
Such correlations persist for RL blue outliers. In all four 

22 V([0 III]) are 220, 400, 480 and 520 kms" 1 for J1634+4809, 
J1505+0326, J1305+5116, and J1443+4725, respectively 
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FIG. 12. — Ratio of the Fe II multiplets to the total H/3 line flux, R 4S10 
versus H/3 line-width in FWHM for the radio-loud NLS1 galaxies (fill cir- 
cles), as well as normal NLS 1 galaxies in the Zhou et al. (2006) sample (small 
dots). 



FIG. 14. — Distribution of Fe II/H/3 for the radio-loud NLS 1 galaxies (solid 
line). For comparison, also plotted is the distribution for the normal NLS1 
galaxies in the Zhou et al. (2006) sample, normalized to the size of the radio- 
loud sample (dashed line). 
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FIG . 1 3 . — Ratio of the Fe II multiplets to the total H/3 flux, R4570 versus the 
ratio of the [O III]-to-H/3 line flux. Plot symbols are the same as in Figure 12. 

we detect blueshifted [Ne V] lines with velocity shifts up to 
-2300 kms" 1 (J1443+4725). The high frequency of blue out- 
liers in RL NLS1 galaxies may indicate that jet-cloud interac- 
tions contribute to, or dominate, the [O III] A5007 blueshifts in 
the RL sources. 

Most of the objects show symmetric profiles of the 
[O III]A5007 line. In three objects, particularly J1633+4718, 
the [O III] doublets show a blue wing that can be modeled by 
an extra Gaussian component (see Figure 1). This indicates 
the presence of possible outflows associated with the narrow 
line region. In addition, J103 123.73+423439.3 shows double- 
peaked profiles in both [O III] A5007 and H/3. 

5. HOST GALAXIES 

The host galaxies of radio-loud NLS1 AGN are poorly 
understood due to very sparse observations. The only ob- 



ject, as far as we are aware of, observed with the HST is 
2MASXJ0324+3410. Its HST image reveals a ring or one- 
armed galaxy morphology, indicating a possible galaxy inter- 
acting/merging history (Zhou et al. 2007). As for our current 
sample, given the relatively high redshifts, most of their host 
galaxies cannot be examined with the SDSS images. The ex- 
ceptions are J1633+4718 and J1644+2619, which have red- 
shifts of 0.1 16 and 0.144, respectively. Interestingly, the host 
galaxy of J1633+4718 has a size of more than 20kpc along 
the major axis and is in an interacting pair/merger system, in 
which the second galaxy hosts a reddened starburst nucleus 
(Bade et al. 1995). Both of the member galaxies appear to be 
gas rich with disturbed morphology, and thus of late morpho- 
logical type. The host galaxy of J 1644+26 19, though being 
compact (<~ 15 kpc), also seems to be disk-like. We note that 
the black hole masses in both galaxies are small, 2. Ox 10 6 and 
8.4xlO 6 M for J1633+4718 and J1644+2619, respectively 
(see § 6.4). 

The SDSS spectrum of J1633+4718 includes a significant 
contribution of the host galaxy light within the fiber of 3" 
diameter, which appears to be dominated by a young stel- 
lar population, as can be seen in Figure 1 (the contamina- 
tion from the starburst nucleus of 4" away should be negligi- 
ble). Besides, marginal imprints of high-order Balmer absorp- 
tion lines, though being week, are likely present in the SDSS 
spectra of other three objects at higher redshifts (J 1 146+3236, 
J1505+0326, and J1634+4809; z = 0.4-0.5), indicating pos- 
sible young stellar populations in their host galaxies. 

We thus suggest that some of the radio-loud NLS1 AGN 
(at least those with a small black hole mass; see § 6.4) may 
reside in gas-rich, disk galaxies, and some reside in galaxies 
with possible young stellar populations. This result contrasts 
clearly with classical RLQs and BL Lac objects, that prefer- 
entially reside in giant elliptical or massive bulge-dominated 
galaxies, or strongly passively evolving galaxies (e.g. Bahcall 
et al. 1997; Dunlop et al. 2003; Floyd et al. 2004; O'Dowd 
& Urry 2005). While it is, to some extent, in line with the 
finding of Ho & Peng (2001) and Ho (2002) for the Seyfert 
galaxies in their samples, which are much less luminous and 
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more sub-Eddington than our radio-loud NLS1 AGN stud- 
ied here, however. As is shown in §6.4, most of the ob- 
jects in our sample have relatively small black hole masses 
(M B h^ 10 6-8 M q ). Hence, considering the M B H-spheroid 
mass relation known for nearby galaxies (Magorrian et al. 
1 998), our finding about the host galaxies of radio-loud NLS 1 
AGN is not unexpected. Deeper observations are needed to 
reveal the host galaxy properties for the remaining radio-loud 
NLS 1 AGN with larger M BH , which are not an easy task given 
their relatively high redshifts, however. 

6. INTERPRETATION AND DISCUSSION 

6.1. The sample 

Our objects fulfill the conventional definition of NLS1 
galaxies, and show typical NLS1 characteristics such as 
strong optical Fe II emission. They also share similar dis- 
tributions in black hole mass and Eddington ratio with 'nor- 
mal' NLS1 galaxies (see §6.4). Therefore, they should 
be bona-fide NLS1 galaxies. The 5 GHz radio luminosities 
L 5G hz range from 10 24 17 to 10 26 96 WHz" 1 with a median of 
jq25.23 Thu^ most of our objects would remain qualified 
as radio-loud even if the alternative criterion of RL AGN is 
applied, which is based solely on radio luminosity, Lsghz > 
10 25 WHz" 1 as proposed by some authors (e.g. Miller et al. 
1990; Kellermann et al. 1994). Of our sample, 19 were spec- 
troscopically targeted in the SDSS as quasar candidates based 
on their optical colors, and 3 were targeted serendipitously; 
only one object was selected solely as a ROSAT/FIRST coun- 
terpart. Hence our very RL NLS 1 sample can be regarded as 
an essentially optically-selected one, similar to other SDSS 
radio-quiet NLS 1 samples. 

The sample spans a luminosity range from Seyfert galaxies 
to quasars, with roughly half of the sample in each category. 
Compared to the 'normal' NLS 1 galaxies in the Zhou' 06 sam- 
ple, which have an average H/3 luminosity (Lh/3) = 3.0 x 
10 42 ergs _1 , our RL objects have relatively higher luminosi- 
ties ({Lap) = 6.0 x 10 42 ergs _1 , see also Figure4). This dif- 
ference is likely a consequence of the trend that RL objects 
are more likely to be found in AGN with more massive black 
holes (and thus higher luminosities), as suggested recently 
(e.g. Best et al. 2005). We find that such very RL objects 
represent only 0.7% of NLS1 galaxies drawn from the SDSS, 
though this fraction may be somewhat uncertain given all 
kinds of selection effects introduced in the surveys. 

Compared to the other two previously published RL NLS 1 
samples of Komossa et al. (2006b) and Whalen et al. (2006), 
which comprise 11 and 16 objects respectively, our current 
sample differs mainly in that it has a much higher radio- 
loudness cutoff. Consequently, the sample is biased for select- 
ing fiat-spectrum sources dominated by beamed radio emis- 
sion. In addition, our sample has a statistically higher red- 
shift distribution than those two samples. In comparison, 
only three and four objects from those two samples satisfy 
Ria> 100, respectively, among which, two from the former 
and one from the latter sample are included in our current 
sample. 

6.2. Broad band continuum emission 

6.2.1. Radio emission: relativistic beaming and jets 

The radio luminosities of our sample objects are well above 
those expected from the most radio-luminous starburst galax- 
ies - 10 22 3-23 4 WHz" 1 (Smith et al. 1998). In addition, no 
imprints of vigorous starbursts are found in the optical SDSS 



spectra of sources of our sample; these spectra sample the cen- 
tral region of the host galaxies within a radius of 3-1 1 kpc. 
None of the objects are found to be associated with IRAS 
sources except for J1633+4718, for which the infrared emis- 
sion is expected (at least predominantly) from the starburst 
nucleus 4" away, rather than from the active nucleus, however. 
These results suggest that for essentially all of the objects the 
radio emission is of AGN origin, rather than of starburst ori- 
gin. 

Large amplitude flux variations up to ~ 70% have been ob- 
served in some of our objects with flat radio spectra. To make 
sure that this is not spurious due to contamination of nearby 
sources, we visually examined the FIRST and NVSS images 
of the variable objects. We find that for 6 of the 8 objects 
no other sources are detected within 3'; for the remaining 2, 
faint sources are present but are all well separated from them. 
Therefore, the observed variations are not caused by flux con- 
tamination from nearby sources. In addition, we argue in this 
paper (see Appendix D for details) that those observed varia- 
tions are largely intrinsic, rather than extrinsic that is caused 
by refractive interstellar scintillation (RISS). 

We have shown that most of the sources known to be fiat- 
spectrum in our sample have very high brightness tempera- 
tures 10 11-13 K; in 4 of the objects T B even exceeds the inverse 
Compton limit of 10 12 K. Such high brightness temperatures 
indicate that the radio emission is highly Doppler beamed and 
hence must originate from (at least mildly) relativistic jets, as 
in blazars. Among the 4 with T B > 10 12 K, J0849+5108 is 
a known OVV/HPQ blazar (e.g. Arp et al. 1979; Sitko et al. 
1984, and see Appendix B) and the high brightness temper- 
ature is well expected; in turn, this gives us confidence in 
the results for the other objects. J0948+0022 has been pre- 
viously suggested to be so by Zhou et al. (2003). For the 
remaining two, J 1505+0326 and Jl 548+35 11, the evidence 
for relativistic beaming is presented for the first time. For 
J 1505+0326, this result is consistent with its VLBI images 
made at 2.3 GHz, which are dominated by a very compact 
component, marginally resolved at a few tens milli-arcsec 
resolution (Dallacasa et al. 1998). Besides, two objects, 
J1633+4718, and J1644+2619 have also been suggested to 
be blazar-like based on VLBI observations (Doi et al. 2007). 
Thus, in a total of 6 objects of our sample there is evidence for 
relativistic beaming. Assuming that the brightness tempera- 
tures of 10 11_12 K are probably also an indication of beaming, 
the number of such objects will increase to 9, i.e. all the 8 
flat-spectrum sources and one with unknown radio index. 

Our suggestion for the presence of (at least mildly) rela- 
tivistic jets is also supported by polarization detected in a few 
individual objects, e.g. J0849+5108 (see AppendixB) and 
2MASXJ0324+3410 (Neumann et al. 1994). Furthermore, 
this is consistent with the observed large amplitude variabil- 
ity in the optical in a few objects monitored, e.g. J0849+5108, 
J0948+0022 (Zhou et al. 2003), since normal NLS1 galaxies 
are found to be much less variable in the optical band (Ai et 
al. 2008). 

These results are most naturally explained in terms of 
beamed emission; and thus we postulate that relativistic jets 
are likely present in most, if not all, of the flat-spectrum, 
very radio-loud NLS1 AGN. However, it should be noted that 
part of the observed properties above may also be explained 
by some other non-blazar, compact sources. For instance, 
flat and inverted spectra are also typical of Gigahertz-Peaked 
Spectrum (GPS) sources, though they generally lack of large 
amplitude flux variations (e.g. O'Dea 1998). Moreover, flat 
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spectra are often seen also in the radio sources of some 'nor- 
mal' Seyfert galaxies, though mostly with much less power 
compared to the objects studied here (e.g. Ho & Ulvestad 
2001). Incorporating the properties of these flat-spectrum 
NLS1 AGN in other wave-bands, we regard the blazar-like 
scenario as the most likely origin of their multi-waveband 
(non-thermal) emission in this paper, though other possibili- 
ties cannot be ruled out completely at this stage. We note that 
the majority of the sources of Komossa et al. (2006b) sample 
(R\a> 10) are steep-spectrum sources, in which strong beam- 
ing is generally not expected. The difference is largely at- 
tributed to the much higher radio-loudness cutoff of our sam- 
ple, as mentioned above. 

6.2.2. Optical continuum: dust extinction and jet 
contribution? 

The majority of our sample objects (16/23) have optical 
continuum slopes bluer than a\=-\. If the continuum is in- 
trinsically blue in all objects (assuming a\=-\.l) but is red- 
dened by extinction, a slope of a\ =-1 corresponds to an 
amount of extinction E(B-V)~ 0.2, i.e. a correction factor of 
only <~ 2 in the B-band flux. So, even if this was the case, most 
of our objects would still have R\ 4>100, above the radio- 
loudness cutoff of our sample. There are 5 objects at low red- 
shifts (z < 0.38) having both the Ha and H/3 lines measurable, 
including Jl 138+3653, the only one showing a significant de- 
viation from a power-law optical continuum. Jl 138+3653 has 
a broad line Balmer decrement of Ha/H/3=4.43, indicative of 
mild extinction. Adopting the intrinsic Ha/H/3 ratio of 3.37 
[with a dispersion of 0.06 dex in log(Ho/H/3)] for blue radio- 
loud quasars (free from dust extinction) found in our recent 
work (Dong et al. 2008) and the Galactic extinction curve, 
we find an internal extinction correction As — 0.97 [E(B- 
V)=0.236], and an extinction corrected R\ 4 a factor of 2.4 
lower, leading to an intrinsic radio-loudness R\a— 86. For 
the remaining 4 objects, the broad line Ha/H/3 ratios range 
from 2.8-3.8 with a mean of 3.4, that are well consistent with 
the intrinsic Ha/H/3 distribution and hence with little or no 
extinction; therefore their radio-loudness R\ 4 measurements 
as above are almost un-affected. Particularly, J 103 1+4234 
shows a red optical spectrum with a\ =0.55, indicating that 
in at least some of the objects the red optical spectrum is in- 
trinsic, rather than resulting from dust extinction. Another 
example is the blazar J0849+5108 with the reddest slope, 
aA=-0.20, that is close to the previous measurement at its 
high light state. Its slope steepens with decreasing luminos- 
ity and changes rapidly on timescales on several months, and 
hence the red slope is preferentially explained as being intrin- 
sic 23 (Arp et al. 1979), similar to red BL Lac objects (Stein et 
al. 1976). 

Observations have revealed that in a large fraction of FS- 
RQs the optical light can be contributed significantly by 
the non-thermal emission from relativistic jets that is highly 
beamed (e.g. Whiting et al. 2001). We have discovered that 
the observed optical fluxes are systematically in excess of 
those predicted from the broad H/3 luminosity, presumably 
from an accretion disk, assuming that radio-loud NLS 1 AGN 
have the same conditions in the broad line region as normal 
NLS1 AGN (see §3.2). In addition, the optical slopes of 
our radio-loud objects are systematically bluer than that of 

23 Alternatively, the red slope is explained as contamination from a fore- 
ground red galaxy on the line of sight (Stickel et al. 1989), though we consider 
this unlikely. 



radio-quiet NLS 1 AGN. These facts point to the presence of a 
second component superposed on the normal AGN-like (ther- 
mal) continuum. In light of the similar radio properties of our 
RL NLS1 galaxies to FSRQs, we attribute the excess opti- 
cal emission to arising from the non-thermal radiation from 
jets, as in FSRQs, i.e. the extension to higher energies of 
the radio-to-infrared synchrotron radiation. This model ex- 
plains naturally the observed difference in the observed-to- 
predicted luminosity ratios between flat- and steep-spectrum 
NLS1 AGN. In flat-spectrum sources, the relativistic beam- 
ing effect is strong — as suggested above — and hence the jet 
component over-shines the thermal optical component; while 
in steep-spectrum sources, the former becomes comparable to 
or weaker than the latter. 

The observed optical continua have a variety of slopes, 
ranging from moderately red to blue. They are reminiscent 
of the optical emission of BL Lac objects, which have both 
red and blue continua (e.g. Stein et al. 1976), corresponding 
to the low- and high-peak frequency of the synchrotron hump 
of LBLs and HBLs, respectively. In fact, this comparison is 
particularly relevant in light of our finding that the RL NLS 1 
AGN have both LBL- and HBL-like SEDs: the blue optical 
continuum may be explained as contribution from the syn- 
chrotron emission peaked at high energies (UV or higher); 
while the red slope may be contributed by that peaked in the 
infrared. In turn, this is consistent with our above inference 
that the red slopes in our objects are mostly intrinsic, rather 
than resulting from extinction. 

6.2.3. X-ray radiation and broad band SED 

So far, we have shown that relativistically beamed jet emis- 
sion, as found in blazars, provides a good explanation of 
the observed radio and optical properties in flat-spectrum RL 
NLS1 AGN. This scenario is also consistent with the fact that 
some of the objects show SEDs similar to FSRQs/SSRQs, 
while the others similar to HFSRQs in the a ro , a ox , and a m pa- 
rameter space (Figure 9). In fact, a few of these objects were 
simply classified as FSRQs or SSRQs in previous studies. It is 
thus natural to speculate that the similar radiation mechanisms 
are operating in these objects as in normal FSRQs/SSRQs and 
HFSRQs for producing the broad-band non-thermal emission. 

From the soft X-ray spectral shapes alone it is not possible 
to distinguish the origin of the X-ray emission, i.e. that asso- 
ciated with accretion disk as in normal RQ NLS 1 galaxies or 
non-thermal jet emission as in blazars, since RQ NLS1 galax- 
ies also show a similar T distribution with a large scatter (e.g. 
Boiler et al. 1996; Grupe et al. 1999; Leighly 2000; Zhou et 
al. 2006). The only exception is perhaps the flat hard X-ray 
component (r = 1 .37 ± 0.49) in J 1633+47 18, whose best esti- 
mate value is hardly found in RQ NLS 1 galaxies but is typical 
of FSRQs. 

One way is to search for possible enhancement in the X- 
ray emission of RL NLS 1 AGN compared to RQ objects in 
a statistical sense, which, if is the case, may be considered 
as an indirect evidence for the presence of a non-thermal jet 
X-ray component. Here, we simply use the X-ray detection 
rate in the RASS as an indicator of possible enhancement in 
the X-ray emission of our objects over normal NLS 1 galaxies 
by comparing the RL and RQ samples. This is plausible since 
our RL NLS 1 sample is essentially an optically selected one 
(see discussion in § 6.1), as the SDSS RQ NLS1 samples. The 
RASS detection rate for RQ NLS1 galaxies is ss 31% from 
the SDSS DR3 NLS1 sample (Zhou et al. 2006); while it is 
50(±1 1)% (1 1/22) for the present RL sample (excluding one 
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solely RASS-selected target). If we consider the objects (both 
detections and upper limits) falling within the HBL locus only 
(Figure 9; excluding 2MASX J0324+34 1 0), the detection rate 
is even higher 58(±14)% (7/12). Therefore, the higher RASS 
detection rate in our RL sample compared to the RQ NLS1 
sample is suggestive of the presence of a significant excess 
X-ray component, presumably from the relativistic jets. 

In Appendix E, we discuss in detail the possible mecha- 
nisms of the X-ray radiation of our RL NLS1 AGN, which 
are summarised here. In all of the objects, the X-rays pro- 
duced in the accretion process as in normal RQ NLS 1 galax- 
ies are likely not negligible and make a considerable contri- 
bution. In addition, we suggest that the non-thermal X-ray 
emission from relativistic jets may also play a similarly im- 
portant role. Such a component is most likely to originate 
from inverse Compton radiation and has a flat X-ray spec- 
trum (in J1633+4718 and possibly in J0849+5108), as seen 
in FSRQs, for objects having the LFSRQ-like SEDs (NLS1 
LFSRQs); whereas it is possibly the high-energy tail of syn- 
chrotron radiation with a steep spectrum for objects with the 
HFSRQ-like SEDs (NLS 1 HFSRQs) . 

6.2.4. Dichotomy ofLFSRQ- and HFSRQ-like NLS1 AGN 

In Figure 9, there are about 7-12 objects, or 1/3-1/2 
of the sample, falling within the box-shaped locus and 
thus having HBL/HFSRQ-like SEDs (the uncertainty comes 
from those with lower a ox limits). Considering that 
radio/optically-selected ('broad' -line) FSRQs have predom- 
inately LBL/LFSRQ-like SED, such a high fraction of 
HFSRQ-like objects is surprising for our sample, which is es- 
sentially optically/radio-selected. The figure also shows that 
there are no objects resembling extreme HBL/HFSRQs (the 
lower-left corner of the HBL locus), which have the syn- 
chrotron component peaked in the X-ray regime. In fact, 
nearly all the candidate HFSRQ-type SEDs have the syn- 
chrotron humps peaked likely around the UV band. 

Strong emission lined-FSRQs with HBL-type SEDs have 
recently been reported as a result of the identification 
of radio-X-ray (ROSAT) selected sources (Padovani et al. 
1997b; Perlman et al. 1998; Padovani et al. 2003). One 
major distinction with our objects is that those objects have 
much broader lines with FWHM > 2000 km s" 1 , typical of 
classical FSRQs. The fraction of HBLs/HFSRQs is sample- 
selection dependent; the authors proposed a fraction (HBLs + 
HFSRQs) of about 10-15% among the whole blazar popula- 
tion and about 30% in the RASS-Green Bank (RGB) sample 
(Padovani et al. 2003). The RGB sample is comparable in 
the way of selection to our RASS-detected sub-sample here, 
which comprises 12 objects; among them, 8 objects are pos- 
sibly of the HFSRQ-type (Figure 9), making the fraction as 
high as ~ 2/3. It seems that the HFSRQ fraction of our radio- 
loud NLS 1 sample is higher than that of HFSRQs in Padovani 
et al. (2003). 

Within the unified blazar sequence, the energy of the syn- 
chrotron peak was found to be well anti-correlated with 5 GHz 
radio luminosity 24 , and with a rx , which was suggested to be 
a tracer of the peak energy (Fossati et al. 1998). We plot in 
Figure 15 a IX versus 5 GHz radio luminosity for the sample, 
with those objects without X-ray detection indicated by ar- 
rows. It clearly reveals a good correlation; the Spearman cor- 
relation test incorporating lower limits of a ra using ASURV 
(Isobe et al. 1986; Isobe & Feigelson 1990) yields a signif- 

24 ^5GHz is thought to be a good tracer of the bolometric luminosity. 
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FIG. 15. — Relationship of a rx and 5 GHz radio luminosity. Symbols: ar- 
rows indicate a rx upper limits for objects without X-ray detection; filled dots: 
flat spectrum radio sources; open circles: steep spectrum radio sources; aster- 
isks: sources without available radio spectral indices. 2MASX J0324+3410 
is also plotted (square). The dotted line corresponds to the simple division 
between LBLs and HBLs using a rx . 

icant correlation at the probability level Pch mce =0A%. ^ can 
be seen that the transition between LFSRQs and HFSRQs (us- 
ing a rx « 0.78) happens at around vLv <~ 10 42 ergs _1 , con- 
sistent with what was previously found between HBLs and 
LBLs (Fossati et al. 1998, see their Figure 7). If a TX does 
trace the synchrotron peak energy, this result shows that our 
radio-loud NLS1 AGN do follow the blazar sequence well. 
In fact, the objects at the low luminosity end are among the 
least radio luminous blazars known so far, falling close to or 
even below the low-end of the FSRQ luminosity function ac- 
cording to unification schemes (Lsghz ~ 10 31,5 ergs" 1 Hz" 1 , 
Urry & Padovani 1995). However, there seems to be no or 
at most a weak correlation between a TX and H/3 luminosity 
(P c s £ ance =0.1). Such a blazar sequence also explains naturally 
the high fraction of HFSRQs, as the majority of our sample 
have lower radio powers compared to classical FSRQs. In 
fact, the objects of our sample have intrinsic powers mostly in 
the range for FR I type radio sources (see § 6.3). 

6.3. Intrinsic radio power and radio-loudness 

Optically thin radio emission at low frequencies is thought 
to be isotropic, and thus can be used as a measure of the in- 
trinsic power of the radio sources free from Doppler boost- 
ing. The low frequency power of radio galaxies can be sepa- 
rated into two distinct classes at the fiducial value Lhsmhz — 
2 x 10 25 WHz" 1 sr" 1 : the FRIs with the radio emission de- 
creasing outwards away from the nucleus, and the FRIIs 
which have lobes with prominent hot spots and bright edges 
(Fanaroff & Riley 1974). It has been suggested that RLQs 
are found in FR II sources, and the distinction of the FR I and 
FR II types is likely determined by the accretion process of 
the central engine (e.g. Baum et al. 1995), as well as by the 
density of the environment (Gopal-Krishna & Wiita 2002). 
Despite a lack of adequate observations on their radio mor- 
phology, it would be informative to compare the isotropic ra- 
dio power of our RL NLS1 galaxies with the FR types. 

We estimate the rest frame 151 MHz luminosities, Lisimhz, 
for our objects, based on their radio data at low frequencies, 
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which are summarised in Appendix C. For 4 objects with 
151 MHz fluxes, Lisimhz are calculated directly by using the 
151-327MHz indices for K-correction. Among the rest, 11 
objects have fluxes at 327 MHz, for which Lisimhz are esti- 
mated by extrapolating the spectra down to 151 MHz assum- 
ing the average aj5™=-0.34 found from the above 4. Using 
the steepest (-0.62) and the flattest (+0.09) a^™ found above 
would increase and decrease Lisimhz by a factor of up to 1.4 
and 1.8, respectively; and thus the estimated Lisimhz are rela- 
tively insensitive to the choice of a 15]^. For the 15 objects 
with either detected or estimated 151 MHz fluxes, Lisimhz 
range from 4.5 x 10 24 to 3.5 x 10 27 W Hz" 1 , with a median of 
7.1 x 10 25 WHz -1 . Among these, three have luminosities in 
the FRII range (J1047+4725, J1305+51 16, and J1443+4725, 
all detected in the 151 MHz survey); while the remaining ma- 
jority have similar luminosities as FRIs or are close to the 
transition line. With estimated Lisimhz =9 x 10 24 WHz _1 , 
2MASXJ0324+3410 also has a luminosity typical of a FRI 
galaxy. However, it is not clear whether the clear separation 
in both radio power and morphology for radio galaxies is ap- 
plicable to RL NLS1 galaxies, considering their apparently 
peculiar, compact morphology. 

It has been suggested that, in classical RL AGN, isotropic 
radio emission at low frequencies can be used as a good in- 
dicator of the power of the jets (e.g. Rawlings & Saunders 
1991). If this is applicable to RL NLS1 galaxies, the above 
results imply that our sample has much lower jet power com- 
pared to classical RLQs, which have radio power and mor- 
phology predominantly in the FRII category (e.g. Rawlings 
& Saunders 1991; Baum et al. 1995; Liu et al. 2006). We will 
address this issue in a detailed study in the future. 

The radio-loudness, by its definition, is largely an apparent 
rather than intrinsic parameter for the sources in which the 
beaming effect is important above 1 GHz in the radio but less 
so in the optical band. As a consequence, some radio-loud 
AGN would have been classified as radio-intermediate or even 
radio-quiet, which might also possess weak and mildly rela- 
tivistic jets in at least some (e.g. Falcke et al. 2006b; Blundell 
et al. 2003; Brunthaler et al. 2005; Wang et al. 2006), were 
their jets not directed close to our line-of-sight (Falcke et al. 
2006a; Jarvis & Mclure 2002). A more physically meaning- 
ful parameter is the 'intrinsic' radio-loudness free from the 
beaming effect. We calculate the intrinsic radio-loudness us- 
ing the un-beamed radio component, that can be extrapolated 
from a low-frequency flux and dominates fluxes at several 
GHz after the jet emission is de-boosted. The intrinsic radio- 
loudness at 5 GHz, ^5c r Hz , is estimated by using 327 MHz 
fluxes for those having this flux measurement. We find R 1 "q Hz 
in the range 34-2600 (with a median of 1 10) assuming a spec- 
tral slope of -0.7 for the un-beamed component, or 17-1380 
(with a median of 56) assuming a slope of -1. It should be 
noted that such estimates are subject to large uncertainties, 
partly due to the above finding that the optical light is also 
partially beamed in some of the objects. Nevertheless, they 
seem to indicate that a sizable fraction of the sample, includ- 
ing 2MASXJ0324+3410, are likely radio-intermediate (e.g. 
<250 was suggested by Falcke et al. 2006a, for flat-spectrum 
sources). Among the sample, there are almost no genuine in- 
trinsically RQ objects with R <C 10, however. 

6.4. Black hole masses of radio-loud NLS1 galaxies 
6.4.1. Estimation of black hole mass and Eddington ratio 



It has become possible recently to estimate the black hole 
masses Mbh of AGN from single epoch spectroscopic obser- 
vational data (e.g. Kaspi et al. 2000; Wandel et al. 1999), 
though the uncertainties of both systematics and scatter re- 
main large (Vestergaard & Peterson 2006; Collin et al. 2006). 
Assuming that the motion of the line-emitting clouds are viri- 
alised, the black hole mass can be determined by Mbh = 
fR BLR Av 2 /G, where / is a scale factor depending on the 
kinematics and geometry of the BLR, Av the emission line 
width, and Rblr the BLR radius. We assume that the com- 
monly used /?BLR-luminosity relation, derived via reverbera- 
tion mapping for nearby, normal AGN, is applicable to NLS 1 
AGN. We use the relation for the H/3 luminosity given by 
Kaspi et al. (2005, their Eq. 1 and Table 3); we intend not to 
use the continuum 5 100 A luminosity for possible contami- 
nation from jet emission (see Section 3.2). A common prac- 
tice, among a few others 25 , is to take FWHM as Av and 
/ = (V3/2) 2 = 0.75, assuming an isotropic distribution with 
random orbital inclinations of the broad line clouds (Netzer 
1990). We follow the practice in this work, the same as in 
Komossa et al. (2006b). 

The estimated black hole masses of our sample range from 

1 x 10 6 to 1.5 x 10 8 M Q , peaked at near the median 1.3 x 
10 7 M Q . The black hole mass distribution is well consis- 
tent with that of the Komossa et al. (2006b) sample. How- 
ever, by studying the systematic effects in Mbh estimation, 
Collin et al. (2006) argued that the use of FWHM rather 
than line-dispersion (crime) would introduce systematic bias in 
the scaling factor /. We also give alternative M B h estima- 
tion by correcting this effect using the formalism proposed by 
Collin et al. (2006, their eq.7). For the objects in our sam- 
ple (FWHM< 2200 km s" 1 ) specifically, this correction sets a 
value of the scale factor / ~1.5, about twice the value in the 
former estimation. This results in M B h estimates a factor of 

2 larger than the values obtained above, ranging from 2 x 10 6 
to 3 x 10 8 M Q with a median 2.6 x 1O 7 M . The black hole 
masses thus estimated are listed in Table 4, which are used in 
the following analysis. 

Having estimated the masses of the central black holes, now 
we can derive the Eddington ratios, Rem- For bolometric lu- 
minosity correction, we assume Lb i=9ALA5ioo, as suggested 
by, e.g. Warner et al. (2004), Vestergaard (2004), and Kaspi 
et al. (2000). The ALasioo luminosity is estimated from the 
broad H/3 line luminosity using Eq.(5) of (Zhou et al. 2006) to 
eliminate the effect of the extra jet contribution. It should be 
noted that such an estimate is subject to large systematic un- 
certainties inherent in estimation of Mbh, as well as of bolo- 
metric luminosity. Using the M B h values obtained with the 
scale factor / = 0.75 (see above), we find the Eddington ratios 
ranging from /? E dd =0.52 to 2.6 and peaked around a mean of 
1.0, similar to those found by Komossa et al. (2006b) with a 
mean of 1.2. However, when the correction of Collin et al. 
(2006) is applied, the Mbh values become about twice larger 
and, consequently, the obtained /? E dd values become about 
twice smaller, i.e. ranging from 0.27 to 1.28 with an average 
of 0.52. 

Now we compare the distributions of Mbh and Rem be- 
tween our RL sample and the normal NLS1 sample of 
Zhou'06. Our RL sample has significantly higher M B h and 

25 A few other scaling relations for Mbh estimation were recently pre- 
sented in the literature(e.g. Onken et al. 2004; Greene & Ho 2005, 2007; 
Vestergaard & Peterson 2006; Collin et al. 2006), that slightly differ from 
one another. 
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/?Edd distributions compared to the RQ one. This agrees with 
the general trend that, as aforementioned in § 6. 1, RL objects 
are more likely to be found in AGN with massive black holes 
(e.g. Best et al. 2005). Given their narrow ranges of the line- 
width, the black hole masses of NLS 1 samples, as determined 
from the H/3/Ha line-width and luminosity, is strongly corre- 
lated with luminosity, and in turn, with redshift for flux lim- 
ited samples. This may explain the higher redshift distribution 
of our RL sample compared to that of the normal NLS 1 sam- 
ple of Zhou'06 (see Figure 2). 

It should be pointed out that the virial M B h estimation is 
subject to assumption of an isotropic geometry of the BLR 
with random orbital inclination of the clouds. If the BLR 
is flat, a small line width would be a result of small incli- 
nation, and then M B h would be largely underestimated (see, 
e.g. Wills & Browne 1986; Decarli et al. 2008). Indeed there 
are some suggestions in the literature for a flattened geometry 
of the BLR in some RL AGN (e.g. Vestergaard et al. 2000; 
Sulentic et al. 2003; Jarvis & McLure 2006); but see Punsly 
(2007). This caveat concerns particularly those of our ob- 
jects, that have compact radio morphology and those that are 
blazar-like, indicating a small viewing angle. Given that the 
Mbh distribution of our sample is indistinguishable from that 
of the RQ NLS 1 sample (within the same redshift range), this 
comes back to the old question whether NLS1 galaxies, as 
a population, have flattened BLRs and are viewed preferen- 
tially face-on, rather than having small Mbh (e-g- Osterbrock 
& Pogge 1985; Puchnarewicz et al. 1992; Bian & Zhao 2004); 
but see e.g. Smith et al. (2004). 

If RL NLS1 galaxies do possess relativistic jets and hence 
anisotropic radio emission, as we have shown here, they can 
be used as a probe to this question since their inclination 
can be inferred from radio observations in a statistic manner. 
However, with the relatively few data available currently, we 
do not find any statistically significant trends of the depen- 
dence of line-width (Mbh) on probable inclination indicators, 
such as radio-loudness or radio slope. Therefore, though there 
is no compelling evidence supporting the planar BLR nature 
of RL NLS 1 from the present data, we cannot rule out this 
possibility, either. A much larger RL NLS1 sample covering 
a wide range of radio-loudness is desirable to shed light to this 
important question in the future. 

6.4.2. Radio-loudness versus black hole mass 

As one of the fundamental relations in understanding AGN 
jets, the dependence of radio-loudness R on M B h has been 
studied extensively in recent years, with the advent of the easy 
estimation of Mbh in AGN. A R-Mbh correlation was pro- 
posed by Laor (2000) that all RLQs in the PG quasar sample 
have M B h > 3 x 10 8 M q , while the radio-quiet quasars have 
a much lower M B h distribution. This trend was confirmed in 
some studies (e.g. Lacy 2001; Mclure & Jarvis 2004; Metcalf 
& Magliocchetti 2006), but was found to be weak or absent in 
others (Oshlack et al. 2002; Woo & Urry 2002). Furthermore, 
by including Seyferts and low luminosity AGN, Ho (2002) 
found that R is mostly correlated with R EM rather than M B h- 
In a recent comprehensive treatment of this question, Sikora 
et al. (2007) demonstrated that R depends on both R E dd and 
Mbh, with two distinct inverse R-R E dd dependence sequences 
for AGN hosted by elliptical (systematically higher R) and 
disk galaxies (systematically lower R), respectively. Regard- 
less whether there is indeed a trend of dependence or not, most 
aforementioned studies show that, on the M B h-^ plane, there 
is a scarcity of RL AGN with low M B h- Even for those low- 
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FIG. 16. — Dependence of radio-loudness on black hole mass for the radio- 
loud NLS1 galaxies of our sample (filled dots) and those from the sample 
of Komossa et al. (2006b, open squares). Over-plotted are schematic areas 
populated by the bulk of the radio sources of the samples of Laor (2000, solid 
line) and Lacy (2001, dashed line). We note that a small number of (normal) 
RL AGN with small Mbh and large R has also been claimed by Woo & Urry 
(2002), which lies mostly at higher R ~ 10 3 ~ 4 above ours in the diagram; 
however, there remains a controversy as to the estimation of those Mbh (see, 
e.g. Jarvis & Mclure 2002; Laor 2003). 

Mbh RL AGN suggested in the literature (e.g. Oshlack et al. 
2002; Woo & Urry 2002), there remains a controversy as to 
the estimation of M B h and R (see, e.g. Jarvis & Mclure 2002; 
Laor 2003). 

Obviously, as first pointed out by Wang et al. (2001) and 
Zhou et al. (2003), radio-loud NLS1 AGN do not follow the 
above trend found for normal RL AGN (provided the virial 
Mbh as estimated above are largely correct). Moreover, Ko- 
mossa et al. (2006b) demonstrated explicitly that RL NLS1 
galaxies are located in a region of the M B h-^ diagram that 
is sparsely populated by other types of AGN (R ~ 10 1 " 3 and 
M BH <1O 8 M ). 

Similarly, we plot #1.4 versus M B h for our very RL NLS1 
galaxies in Figure 16. Also plotted are the loci of the ob- 
jects in the samples of Laor (2000) and Lacy (2001) for com- 
parison. Two pieces of results can be inferred immediately 
from Figure 16. Firstly, within our sample, there is no cor- 
relation found between Ri 4 and Mbh, nor between Ri 4 and 
•^Edd (^ > chance = 0-44 and 0.27, respectively); however, we note 
that these results are rather uncertain given the small sam- 
ple size and the limited ranges in both R14 and R^d- Sec- 
ondly, our new result confirms and strengthens that of Ko- 
mossa et al. (2006b), and improves upon it by extending the 
RL NLS 1 locus to even lower M B h and higher R 1.4. In partic- 
ular, the bulk of our RL NLS 1 galaxies lies within a locus of 
M B h — 1O 6_7 - 5 M and R ~ 10 1 " 3 , which appears to be pref- 
erentially avoided by other types of RL AGN [compared to, 
e.g. Sikora et al. (2007), see their Figure 4], except for low 
power Seyfert galaxies and LINERs; however, those excep- 
tions have exclusively low Eddington ratios R E dd < 10~ 2 and 
are drastically different from NLS 1 with high R E dd (see also 
Ho 2002). Compared to Figure 4 in Sikora et al. (2007), our 
RL NLS 1 galaxies appear to have the largest radio-loudness 
up to #1.4 ~ 10 3 " 4 among all types of AGN within this M B h 
range. This means that the RL NLS1 galaxies show strong 
excess radio emission with respect to their M B h, when the 
suggested radio power-M B H relation for normal AGN is com- 
pared 26 (e.g. Lacy 2001). Finally, we note that the above 

26 We note that several (normal) RL AGN with small Mbh similar to ours 
and with even larger R (~ 10 3-4 ) have also been claimed by Woo & Urry 
(2002), though there remains a controversy as to the estimation of Mbh (see, 
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two results remain valid even if the Doppler de-boosted radio- 
loudness values (see § 6.3) are used instead, which are in the 
range of R ~ 10 1-3 . 

Hence, the existence of RL NLS 1 galaxies as a population 
with small M B h, if estimated correctly, fills up the previously 
sparsely populated region in the M Bli -R parameter space for 
normal RL AGN. This result indicates that there indeed exist 
AGN with high radio-loudness (R <~ 10 2-3 ) but intermediate 
Mbh (1O 6_8 M ). Though there is the systematic uncertainty 
in estimating M B h ascribed to the possible BLR inclination ef- 
fect, we argue that this result may not be affected qualitatively. 
This is because the Komossa et al. sample, though composed 
of mainly steep spectrum sources (thus relatively larger incli- 
nation), led to similar results. Furthermore, in a few objects 
where the host galaxies are resolved, their small M B h are sup- 
ported by a lack of massive galactic bulges, as predicted by 
the Mbh -bulge relation (e.g. Magorrian et al. 1998). Our re- 
sult also implies that the inverse relation between R and Rem 
as suggested previously (see, e.g. Greene et al. 2006; Sikora 
et al. 2007) may not work at the high R E dd end close to the Ed- 
dington rate. Furthermore, Sikora et al. (2007) argued that RL 
AGN hosted by elliptical galaxies are much radio-louder than 
AGN hosted by disk galaxies. Though this statement may 
be true in general, it is perhaps more appropriate in terms of 
RL fraction, i.e. the fraction of RL objects is lower in AGN 
hosted by disk galaxies compared to that of AGN hosted by 
ellipticals; in addition, the difference in radio-loudness is now 
less than what was thought previously (Sikora et al. 2007). 

6.5. On the nature of radio-loud NLS 1 galaxies 

6.5.1. A previously unknown population of radio-loud AGN? 

On the black hole mass-radio-loudness plane our RL NLS 1 
galaxies form a natural extension to low M B h of classical 
RLQ/BLRGs (Broad Line Radio Galaxies) with roughly con- 
stant radio-loudness (Figure 16; see also Figure 3 in Sikora et 
al. 2007). This is not surprising, since the transition between 
BLAGN and NLS1 AGN is smooth in the line-width of the 
broad lines, and, as a consequence, in black hole mass and 
luminosity. If the Mbh are underestimated systematically (see 
§ 6.4), the loci of the radio-loud NLS1 AGN will move even 
closer to those of normal RL AGN. In this regard, it is natural 
to treat RL NLS 1 galaxies as the 1ow-M B h end of normal RL 
AGN. 

On the other hand, however, RL NLS 1 galaxies show some 
peculiarities in observed properties that are unusual in com- 
parison with, or at an extreme of, normal RL AGN. The most 
distinctive features are in the radio properties. Even for steep- 
spectrum sources, the radio sources of RL NLS 1 galaxies ap- 
pear to be compact, with the size less than several tens of kpc 
or even smaller wherever higher resolution observations are 
available. In contrast, the radio sources of RLQ/BLRGs are 
about 100-400 kpc in a similar redshift and luminosity range 
as our radio-loud NLS 1 AGN (Singal 1993). Moreover, most 
of the RL NLS1 galaxies have radio power in the range of 
FR I radio galaxies, whereas RLQ/BLRGs tend to be associ- 
ated with FRII type sources 27 (e.g. Laing et al. 1983). 

In the following aspects the two types are not mutually ex- 
clusive but overlap to some extent, yet they seem to differ in 
a systematic way. Firstly, There seems to be a higher fraction 

e.g. Laor 2003). Our sample appears to differ from theirs in the sense that 
it fills up the gap between R ~ 10 and 10 3 in the small Mbh regime in their 
Mru-R diagram. 

27 There are exceptions, though; see, e.g. Blundell & Rawlings (2001). 



of HBL/HFSRQ-like blazars in optically selected RL NLS1 
galaxies than in normal RLQ/BLRGs of radio/optically se- 
lected samples (Padovani et al. 2003), though further con- 
firmation using complete samples is needed. Secondly, as 
commonly suggested, NLS1 AGN are accreting at close to, 
or larger than the Eddington limit (e.g. Collin & Kawaguchi 

2004) . At such high rates the accretion flow proceeds via a 
slim disk where energy is partially advected into the black 
hole before being radiated due to relatively long cooling time 
(e.g. Abramowicz et al. 1988; Mineshige et al. 2000; Wang 
& Netzer 2003). Thus the accretion processes in radio- 
loud NLS 1 AGN may differ drastically from those in normal 
RLQ/BLRGs, where a standard thin disk is believed to be op- 
erating, though they do overlap at /?Edd — 1 to some extent 
(see e.g. Sikora et al. 2007). The final point — on the host 
galaxies — is only suggestive because of a lack of sufficient 
evidence. RL NLS1 nuclei, at least some, seem to be hosted 
in disk or interacting galaxies, and/or galaxies with relatively 
young stellar populations. In contrast, classical RLQ/BLRGs 
28 are preferentially hosted in giant elliptical/massive bulge 
dominated galaxies or at least strongly passively evolving 
galaxies (e.g. Bahcall et al. 1997; Dunlop et al. 2003; Floyd et 
al. 2004; O'Dowd & Urry 2005), although there was also 
evidence reported for star-formation in some hosts beyond 
z > 0.4 (e.g. Orndahl & Ronnback 2005; O'Dowd & Urry 

2005) . 

In consideration of their unusual properties, we are left with 
two options regarding the nature of RL NLS1 galaxies: ei- 
ther they are a downsized extension of normal RL AGN to 
the lower-M B H (and perhaps higher-/? E dd) range, or they form 
a previously unknown population distinct from normal RL 
AGN. The former would invoke some unknown processes that 
can explain the observed differences between RL NLS 1 and 
normal RL AGN with the smooth change of Mbh- For the 
latter, the operating physical processes in likely black hole 
accretion and jet formation may be radically different for the 
two classes. A related, interesting question is whether the 
differences in the observed properties between RL NLS 1 and 
normal RL AGN can be naturally explained by the same un- 
derlying physical processes that make NLS1 different from 
normal AGN, or whether some additional mechanisms are 
needed. To further understand these questions, a large sam- 
ple of RL NLS 1 galaxies with the full range of radio-loudness 
is needed, as well as extensive observations, especially in the 
radio band. As complementary discussion, we also compare 
them with two special types of RL AGN below, namely, BL 
Lac objects and compact steep-spectrum sources. As will be 
seen, RL NLS 1 galaxies may be a heterogeneous population 
themselves. 

6.5.2. Relationship with BL Lac objects 

One object in our sample, J0849+5108 (0846+5 1W1), was 
formerly classified as a BL Lac object at its high light state 
based on previous observations (see Appendix B for details). 
As such, it presents an interesting case linking possibly RL 
NLS1 galaxies with BL Lac objects, as pointed out by Zhou 
et al. (2005), though its black hole mass (2.6 x 1O 7 M ) is 
lower compared to those typical of BL Lac objects (Mbh ~ 
10 8_9 M Q , e.g. Cao 2003). In fact, it has been speculated in 

28 Ho & Peng (2001) showed that some Seyfert galaxies hosted in disk- 
dominated spirals are in fact RL when the optical nuclear luminosities are 
considered properly. However, the vast majority of their sample objects are 
of substantially low power and sub-Eddington (Ho 2002). 
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the past decade that concealed Seyfert nuclei exist in a few BL 
Lac objects, including BL Lac itself (Vermeulen et al. 1995; 
Corbett et al. 2000). It is therefore tempting to postulate that 
some RL NLS 1 galaxies may actually be BL Lac objects (with 
low M B h) at their low light states. 

The bulk of our RL NLS 1 galaxies have a BLR luminos- 
ity 29 distribution lying within log (Leui/erg s -1 )=43.5-44.5, 
the low end of which overlaps partially with those of BL Lac 
objects (a median of 10 43,7 erg s" 1 was given for the upper lim- 
its of Lblr by Celotti et al. 1997). Specifically, J0849+5108 
has L B lr ~ 10 43 7 ergs" 1 . However, this is likely not the case 
for most of the RL NLS 1 galaxies given their relatively higher 
Lblr compared to the upper limits for BL Lac objects. On 
the other hand, this possibility means that BL Lac objects 
might not be a homogeneous population, and some may just 
be Seyfert nuclei/NLS 1 at a high state of the non-thermal ra- 
diation from the jet. Similarly, this concerns only BL Lac ob- 
jects at the Iow-Mbh end, as is the case of J0849+5108, rather 
than the bulk of classical BL Lac objects that habour massive 
black holes and accrete likely via a low efficiency accretion 
flow (e.g. Cao 2003). 

6.5.3. Relationship with CSS sources 

Compact steep-spectrum sources are compact, powerful ra- 
dio sources with spectral peaks near 100 MHz. They are in- 
herently compact with linear sizes < 20kpc, having (often 
symmetric) double/triple morphology with complete jets and 
lobes on small scales. The most popular scenarios to explain 
their nature advocate that they are either young and evolv- 
ing radio sources or confined by dense gas of their environ- 
ment (see O'Dea 1998, for a review). Very recently, a few 
lines of observational evidence point to an interesting link be- 
tween RL NLS1 galaxies and CSS, which share some degree 
of similarity as the following. (1) As mentioned in § 1, the 
radio sources of all RL NLS1 galaxies are compact, includ- 
ing steep-spectrum sources. In particular, a few NLS1 AGN 
having higher resolution observations show small-size struc- 
ture: e.g. most of the radio emission of 2MASX J0324+34 1 
is from a region smaller than 6.8 mas (8.2 pc projected size), 
as seen from VLB A images at 8.4 GHz (in preparation). (2) 
The low frequency spectrum of J 1047+4725, one of the steep- 
spectrum objects in our sample, appears to flatten at around 
a few hundred MHz as indicated by its flat slope of -0.28 
between 327 and 151 MHz; such a spectral break is charac- 
teristic of CSS. In fact, the RL NLS1-CSS connection was 
advocated by Oshlack et al. (2001) and Gallo et al. (2006) 
for the case of PKS 2004-447, although it is unclear as to 
whether PKS 2004-447 is a genuine NSL1 galaxy (see Zhou 
et al. 2003, for a short comment). (3) There are observational 
indications that the host galaxies of a large fraction of CSS 
are interacting/merging and may contain significant amount 
of dense gas (O'Dea 1998), similar to what was found in some 
RLNLS1 galaxies. 

The radio power distribution of our RL NLS 1 galaxies over- 
laps with the lower part of that of CSS (O'Dea 1998). As 
such, if there is an overlap between CSS and RL NLS1 galax- 
ies, they should be mostly the CSS having low radio power, 
and, relatively small Mbh- Since CSS sources also possess 
(small scale) jets, in cases where the jets are directed close to 

29 The luminosity of the broad line region is estimated from the H/3 lu- 
minosity using the method proposed by Celotti et al. (1997) and the quasar 
emission line ratios given by Francis et al. (1991), yielding Lblr = 25.26Lh/3 
(see also Liu et al. 2006). 



our line-of-sight, we expect to see blazar-like objects. Some 
of those could be the flat-spectrum RL NLS 1 galaxies stud- 
ied here, in which relativistic jets are inferred to be present. 
Conclusive remarks have to await future proper and detailed 
comparisons of the two classes, as well as more radio obser- 
vations with higher resolution, however. 

6.5.4. Implication for jet formation in AGN 

Given the high Eddington ratios of NLS 1 galaxies in gen- 
eral, as commonly suggested and also found here, our RL 
NLS1 AGN present a strong case where jets, possibly at least 
mildly relativistic, can be produced in high R E <id systems. This 
strengthens the similar situation in powerful RLQs with much 
larger M Bll - This contrasts clearly with black hole X-ray bi- 
naries where jets are found at the low/hard state, despite the 
suggestion that transient jets (or mass ejection) may occur at 
the so called very-high state where R Et id is higher (Maccarone 
et al. 2003; Fender 2004). The formation of jets is thought to 
be related to the black hole and accretion process (e.g. Bland- 
ford 2000; Celott & Blandford 2001; Meier 2003; Ballantyne 
2007). The finding of relativistic jets in black hole accre- 
tion systems with high R Edd (where the accretion proceeds 
via probably a new form other than the standard disk) and, 
moreover, with relatively low Mbh is important for address- 
ing basic questions as to how jets are formed. Jet formation in 
high R E dd systems (slim disks) is far from clear, though there 
are some suggestions that outflows as well as mildly beamed 
and collimated emission are likely present in such systems 
(Ohsugaetal.2005). 

As has been shown above (see also Komossa et al. 2006b; 
Whalen et al. 2006), apart from the continuum radiation, 
the RL NLS1 galaxies resemble closely their RQ counter- 
parts in many ways, e.g. optical emission properties (except 
for stronger Fe II emission), the distributions of black hole 
masses and the Eddington ratios. These facts suggest that 
RL NLS 1 galaxies are simply otherwise normal NLS 1 galax- 
ies but having relativistic jets, reminiscent of the RL/RQ di- 
chotomy of more powerful, classical AGN. It would be inter- 
esting to find out whether there is a bimodal distribution in 
the radio-loudness, as reported for the broad-line quasar pop- 
ulation (though this is still a matter of debate), in future work 
considering objects with the whole range of radio-loudness. 
The nearly indistinguishable properties other than the pres- 
ence of radio jets possessed by the two subsets of NLS 1 galax- 
ies, RL and RQ, seemingly hint at a possibly important role of 
black hole spin in the making of relativistic jets in RL NLS1 
galaxies. 

7. SUMMARY OF CONCLUSIONS AND IMPLICATIONS 

A homogeneous sample of genuine RL NLS1 galaxies 
(radio-loudness R\A^t 100) is presented for the first time, 
comprising 23 objects drawn from SDSS spectroscopic ob- 
jects with FIRST detections. Their NLS1 nature is discovered 
by analysing systematically the spectroscopic data from the 
QSO and galaxy database in the SDSS DR5. The sample can 
be treated as an optically-radio selected one. We carry out 
systematic investigations into their multiwaveband properties 
from radio to X-rays using data available from the archives. 
The main results are the following. 

All radio sources of the sample are compact, less than a few 
tens of kilo-parsec at least. Several with VLBI observations 
are unresolved or marginally resolved at scales of about tens 
of parsec. Among those with available radio spectral indices, 
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the majority are flat-spectrum sources (a > -0.5 around 1.4— 
5 GHz), and a few even have inverted spectra (a > 0). Radio 
variability is found to be common in most of the flat-spectrum 
sources. We were able to estimate the brightness temperatures 
of the radio sources, that are as high as T B > 10 11 K and even 
exceeding the inverse Compton limit (10 12 K) in four objects. 
The high brightness temperatures indicate the presence of rel- 
ativistic beaming. The majority of RL NLS1 galaxies has 
a low-frequency (isotropic) radio power falling in the range 
of FRI type radio galaxies (Lisimhz < 10 25 WHz" 1 sr" 1 ), and 
the remaining minority in the range of the FR II type. When 
isotropic, rather than Doppler boosted, emission is consid- 
ered, a large fraction of the objects can be considered as in- 
trinsically radio-intermediate. 

The observed optical luminosities are systematically in ex- 
cess of the thermal continua predicted from the H/3 line lu- 
minosities by a factor of up to several. The optical slopes 
of our RL objects, while spanning a wide range from moder- 
ately red to blue, are systematically bluer than that of radio- 
quiet NLS 1 AGN. We interpret this effect as the contribution 
from a (non-thermal) component, presumably from relativis- 
tic jets, to the optical band superimposing the central ionising 
continuum seen in radio-quiet NLS 1 AGN. The optical Fe II 
emission in the radio-loud NLS1 in our sample is on aver- 
age stronger than 'normal' NLS1 galaxies. Except this, the 
two subsets of NLS 1 galaxies seem to be indistinguishable in 
other emission line properties. 

On the radio-optical-X-ray effective spectral indices plane 
(a ox vs. a ro ), our objects have broad-band SEDs well consis- 
tent with those of blazars, i.e. both the HFSRQ and LFSRQ 
types as well as SSRQs. Of particular interest, the fraction of 
the candidate HFSRQ-type objects seems to be high, though 
their synchrotron component peaks are only moderately high 
(at around the UV). This finding extends the recently discov- 
ered strongly-lined (FSRQ) HBL-type blazars (Padovani et al. 
2003) to NLS 1 galaxies. 

The continuum shape in the soft X-ray band exhibits a large 
variety among our RL NLS1 galaxies. There is evidence in 
one, or possibly two objects for a flat X-ray spectral com- 
ponent (r < 1.8), which we interpret as an inverse-Compton 
origin from jets as in FSRQs. Most of the objects with esti- 
mation of the X-ray spectral slopes show steep spectra. The 
steep X-ray spectrum is also consistent with what is expected 
for HFSRQs, that originates from the high energy end of the 
synchrotron component. Based on this circumstantial evi- 
dence, we suggest that such a component may be a signifi- 
cant contributor to the observed X-ray emission, though the 
X-ray component as seen in radio-quiet NLS 1 galaxies may 
also play a role. 

In view of the independent lines of evidence presented here, 
as well as those found in previous studies for several indi- 
viduals (either inclusive or exclusive of our sample), we sug- 
gest that there likely exists a population of RL NLS 1 galaxies 
showing the broad-band continuum emission characteristic of 
blazars, though alternative explanations cannot be ruled out 
completely. In this scenario, beamed non-thermal emission 
from jets, that are at least mildly relativistic, dominates the 
high-frequency radio band, and contributes significantly to the 
optical and possibly the X-ray bands. Our finding extends the 
radio-loud quasar population into a narrow line-width regime 
of the broad emission lines. We postulate that all genuine 
radio-loud NLS 1 galaxies possess small scale, at least mildly 
relativistic jets, reminiscent of classical radio-loud AGN. In 
this scenario, steep-spectrum NLS 1 sources (perhaps includ- 



ing some of those with smaller radio-loudness within R =10- 
100) are simply objects of the same population whose radio 
jets are aligned at a relatively large angle to the line-of-sight. 

The host galaxies of two nearby RL NLS1 nuclei, whose 
images are resolved, have disk morphology, one even in a 
merger, in line with the finding for the previously published 
blazar-NLSl hybrid object 2MASX J0324+34 1 (Zhou et al. 
2007). In several objects imprints of young stellar populations 
are likely to be present in form of high-order Balmer absorp- 
tion lines in the SDSS spectra. 

Based on the commonly used relation of black hole mass 
with broad line width and luminosity, our RL NLS1 galax- 
ies have estimated Mbh in the range of 10 6 - 10 8 M Q , pro- 
vided that this Mbh scaling relation is applicable to NLS 1 and 
the BLR is not largely flattened. The thus estimated Edding- 
ton ratios range from R^d =0.5 to 3, and peak at /?Edd~ 1 
(assuming Lhoi = 9ALasioo)- RL NLS1 galaxies deviate from 
the known trends of radio-loudness versus the Eddington ra- 
tio and black hole mass, and are populated in a region in the 
radio-loudness-M B H diagram where normal AGN are seldom 
found, i.e. moderate M BH and large R (M BH ~ 1O 6_7 - 5 M , 
R ~ 10 1-3 ), confirming and improving upon the previous re- 
sults suggested by Komossa et al. (2006b). As such, RL NLS 1 
show strong excess radio emission with respect to their M B h, 
when compared to the previously suggested radio power-M B H 
relation for normal AGN (e.g. Lacy 2001). Our results imply 
that large radio-loudness (R = 10 1-3 ) can be produced in AGN 
with moderate Mbh down to <~ 10 6 - 10 7 M Q and accreting 
at/close to the Eddington rate; and some of these AGN are 
hosted in late-type galaxies, some even in interacting systems. 
The presence of relativistic jets in such systems may have in- 
teresting implications for understanding jet formation in AGN 
in general. 

In spite of a likeness to normal RL AGN in many aspects, 
RL NLS1 galaxies show some extreme or unusual proper- 
ties, though some of the evidence are only circumstantial. 
These include compact radio morphology, relatively low ra- 
dio power, a high incidence of HBL/HFSRQ-type SEDs, host 
galaxies of probably late-type and/or with young stellar pop- 
ulations, and the common properties characteristic of NLS 1 
galaxies, e.g. smaller black hole masses, higher Eddington 
ratios, and stronger Fe II emission. All of these lead us 
to postulate that either they are an extension of normal RL 
AGN to the lower-MeH range, and all those differences can 
be explained as a consequence of decreasing M B h; or they 
form a previously unknown population distinct from normal 
radio-loud AGN. As one interesting possibility, we speculate 
that low-power CSS might be the parent population of the 
flat-spectrum RL NLS1 galaxies with beamed jet emission. 
For the future, a much larger sample and extensive multi- 
waveband observations are essential for understanding the na- 
ture of radio-loud NLS 1 galaxies. 

Finally, if the broad-band radiation mechanisms in fiat- 
spectrum radio-loud NLS1 AGN are indeed the same as in 
blazars, we expect some of these objects to be 7-ray emit- 
ters in the GeV band. This can be tested by observations of 
GLAST in the near future. 
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APPENDIX 

A. OPTICAL SPECTRAL MODELING 
The starlight and nuclear continua are modeled as 

S{\)=A host {Eg°Jy,X) A(X)+A nucleus (E n B U -v eus , A) [bB(X) + c b C h (X) + c a C n (X) + dD(X)] (Al) 

where S(X) is the observed spectrum. A(A) = Ym=i a i IQ(\ <7 *) represents the starlight component modeled by our 6 synthesized 
galaxy templates, which had been built up from the spectral template library of Simple Stellar Populations (SSPs) of Bruzual & 
Chariot (2003) using our new method based on the Ensemble Learning Independent Component Analysis (EL-ICA) algorithm. 
The details of the galaxy templates and their applications are presented in (Lu et al. 2006). A(A) was broadened by convolving 
with a Gaussian of width cr* to match the stellar velocity dispersion of the host galaxy; in this process the velocity difference 
between emission and absorption lines is taken into account. The un-reddened nuclear continuum is assumed to be B(X) = A" 1 7 
as given in Francis (1996). We modeled the optical Fe II emission, both broad and narrow, using the spectral data of the Fe II 
multiplets for I Zw I in the AA3535-7530A range provided by Veron-Cetty et al. (2004, Table A1,A2). We assume that the broad 
Fe II lines (Q, in Eq.l) have the same profile as the broad H/3 line, and the narrow Fe II lines (C n ), both permitted and forbidden, 
have the same profile as the that of the narrow H/3 component, or of [O III] A5007 if H/3 is weak. D(X) represents the templates of 
higher-order Balmer emission lines (10 < n < 50) and Balmer continuum generated in the same way as in Dietrich et al. (2003). 
AhnstiEg^yjX) and A nuc i eus (Eg"y eus , X) are the color excesses due to possible extinction of the host galaxy and the nuclear region, 
respectively, assuming the extinction curve for the Small Magellanic Cloud of Pei (1992). The fitting is performed by minimizing 
the x 2 with E^°J V , Eg"y eus , a ; , cr», b, c b , c n , and d being free parameters. 

B. NOTES ON INDIVIDUAL OBJECTS 

J0849+5108 (0846+5 1W1): This object has drawn substantial attention ever since its first identification with a BL Lac object 
(at a bright state) by Arp et al. (1979), at a claimed high redshift of z = 1.86 30 . Its radio flux also showed large amplitude 
variations up to 50% at 1.4 GHz. This object exhibited a remarkable optical brightening with an amplitude of AV <~ 5 within 
less than a month (Arp et al. 1979) and high and variable optical polarization (Moore & Stockman 1981; Sitko et al. 1984), 
making it one of the optically violent variables (OVV) and highly polarized QSOs (HPQ). Night-to-night and even intra-night 
optical variability was also found (Sitko et al. 1984; Stalin et al. 2005). At maximum light outburst, the continuum is featureless, 
whereas its spectra taken at normal/minimum light states show strong emission lines (Arp et al. 1979; Stickel et al. 1989). Thus 
it is a transition between a quasar and a BL Lac object, i.e. appeared to be a BL Lac object only at its flaring state. The large and 
variable polarization in the optical band suggests strongly that this emission is due to synchrotron radiation (Sitko et al. 1984). 
The conspicuous characteristic emission lines identified in the SDSS spectrum of J0849+5108 give a redshift of z = 0.583. A 
detailed study of its SDSS spectrum and other properties was presented by Zhou et al. (2005). The optical spectral slope obtained 
from fitting the SDSS spectrum is red, with a\~ -0.2. 

J0948+0022: The object was first discovered for its radio emission in various radio surveys (Ssghz ^0.2-0.3 Jy). It was 
identified with a NLS1 galaxy, based on its SDSS spectrum, by Williams et al. (2002) and Zhou et al. (2003), and was studied 
in detail as one of the few previously known very RL NLS1 galaxies. The radio source has an inverted radio spectrum and was 
unresolved with VLBI at milli-arcsec resolution. These observations were explained in terms of relativistic beaming (Doi et 
al. 2006). A very high brightness temperature (<~ 10 13 K) inferred from radio flux variations suggest the presence of a jet with 
relativistic beaming (Zhou et al. 2003). Long-term optical variability with an amplitude of about 1 mag in the B-band was also 
noted in Zhou et al. (2003). 

J1505+0326: It has radio fluxes peaked at <~ 7 GHz (Tinti et al. 2005). Its VLBI images are dominated by a very compact 
component, marginally resolved at 2.3 GHz band (Dallacasa et al. 1998). The radio emission at centimeter wavelengths shows 
variability. 

J1633+4718 (RXJ 16333+4719): The object was identified with a NLS1 galaxy in optical identification of ROSAT All-Sky 
Survey (RASS) sources (Moran et al. 1996; Wisotzki & Bade 1997). The AGN (see Figure 1 for its SDSS spectrum) is actually 
one of the two nuclei, separated by 4", residing in an interacting/merger system. The other nucleus has a starburst spectrum 
(Bade et al. 1995). The radio source, which is associated with the AGN, is unresolved with VLBI at milli-arcsec resolution 
(Doi et al. 2007). Its spectral shape above 5 GHz is inverted, as measured by simultaneous multiwavelength observations by 
Neumann et al. (1994). The system was noticed to be associated with the IRAS source IRAS 16319+4725 (RA=16h33m23.03s, 
Dec=+47dl9m03.5s), whose positional uncertainty — much larger than the separation of the two optical nuclei — encompasses 

30 This redshift value, given by the authors based on the spectra with low resolution and signal-to-noise ratios, turned out to be wrong. 
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both of the nuclei (Boiler et al. 1992; Moran et al. 1996; Condon et al. 1998b). The [O III] doublet shows a noticeable blue wing 
in its SDSS spectrum. 

J1644+2619: Its radio source has an inverted spectrum and was unresolved with VLBI at milli-arcsec resolution (Doi et al. 
2007). 

C. RADIO EMISSION AT LOW FREQUENCIES 

In the above sections, evidence has been presented for the blazar-like properties of the (flat-spectrum) RL NLS1 galaxies, 
arguing for the presence of relativistic jets in them. As is known in normal radio quasars/galaxies, radio jets are associated 
with large scale radio lobes. These show diffuse radio emission with steep spectrum (a r ~ -0.77; e.g. Kellermann 1966) which 
dominates the radio fluxes at low frequencies (several hundred MHz), via synchrotron radiation from optically thin plasma. Here 
we examine the radio properties of our sample at low frequencies, though the data are sparse in general, especially at around 
100 MHz. The radio emission is searched for using the WENSS 31 survey (Rengelink et al. 1997) at 327 MHz and using NED 
for surveys at lower frequencies. We find 15 objects with detections at 327 MHz, for which the 1.4 GHz to 327 MHz spectral 
indices can be computed. We found a median a^-m. °f ~0-47, and <~ 47% being steeper than -0.5 with the steepest as 

-0.83, systematically steeper than the 1 .4-5 GHz indices (3/1 1 are steep-spectrum). This means that the spectral indices steepen 
progressively towards low frequencies. 

This behavior can well be explained by onset of a steep-spectrum component dominating the low frequencies, as commonly 
seen in normal RL AGN. Therefore, we postulate that in these objects there exists diffuse, optically thin emission, presumably 
from putative radio lobes that are powered by (relativistic) jets. If this is the case, those blazar-like (flat-spectrum) NLS1 AGN 
must possess the general core-halo morphology as seen in some blazars (e.g. Stannard & Mcilwrath 1982; Browne & Perley 
1986; Murphy et al. 1993). Given the compact morphology unresolved at the FIRST resolution (smaller than several tens kpc), 
there are two possibilities: either they have large scale extended emission — similar to the classical radio sources — but are viewed 
close to the radio axis, or are intrinsically compact, similar to CSS (compact steep-spectrum sources). 

At the even lower frequency of 15 1 MHz, 4 objects have detected radio emission, all from the 6C/7C survey (e.g. Baldwin et al. 
1985; Hales et al. 1988), namely J1047+4725 (1.81Jy), J1305+5116 (0.32Jy), J1443+4725 (0.5 Jy), and J1548+3511 (0.26 Jy). 
The 15 1-327 MHz indices are aj5™=-0.28, -0.54, -0.62, +0.09, respectively, and the average is -0.34. These values are much 
flatter than the 'canonical' slope for diffuse optically thin emission from powerful sources (< -0.7, e.g. Kellermann 1966), 
indicating that the spectra have already turned over, or started to turn over, in this frequency range due possibly to synchrotron 
self-absorption. In fact, the radio spectra of J1047+4725 and J1548+3511 show a clear trend of turnover at low frequencies. 
Considering the relatively large redshifts of the sample, the break frequencies in some sources could be as high as several hundred 
MHz in the source rest frame. Spectral turnovers at such frequencies are often seen in radio sources with small scale structure, 
such as in CSS (e.g. O'Dea & Baum 1997). The previously reported extreme NLSl/blazar hybrid object 2MASX J0324+34 1 
(Zhou et al. 2007) was also detected at 151 MHz (1.02Jy); its radio slope also steepens significantly towards low frequencies, 
a i5iM =- 0-44 (compared to a r ~ 0.1 at high frequencies, Neumann et al. 1994), suggesting the dominance of a steep-spectrum 
component at low frequencies. 

D. EFFECTS OF REFRACTIVE INTERSTELLAR SCINTILLATION ON RADIO FLUX VARIATIONS 

For compact radio sources (milli-arcsec scale), it is known that refractive interstellar scintillation (RISS), caused by inhomo- 
geneities in the interstellar medium, may give rise to variations of radio flux densities (e.g. Rickett 1986, 1990; Blandford et 
al. 1986). RISS is wavelength-dependent and concerns predominantly the low frequency domain (below 1 GHz). At centimeter 
wavelengths, RISS mostly causes "flickering" at a lever of a few percent, as was found in observations (Heeschen 1984) and is 
also consistent with the theoretical models (Rickett 1986; Blandford et al. 1986). Although variations at a larger amplitude were 
sometimes ascribed to RISS, they were limited to a small number of individual objects which ubiquitously are at low Galactic 
latitude (< 10°) or even have a line-of-sight passing through the Galactic plane (Hjellming & Narayan 1986; Dennison et al. 
1987). Monitoring observations of radio sources at frequencies above 1 GHz showed that both intrinsic and extrinsic (RISS) 
variability contribute to observed radio flux variations, with the former having timescales from months to years and the latter 
being faster and having rms amplitude typically < 10% (e.g. Lazio et al. 2001; Mitchell et al. 1994; Rickett et al. 2006). 

In consideration of the large variation amplitudes in some of the objects of our sample and their high Galactic latitudes (b=3 1 — 
51°, Table 3), we argue that the variations are predominantly intrinsic for the majority of the variables reported here, instead 
of arising from RISS. This is especially true for variations at 5 GHz or even higher frequencies. Furthermore, the brightness 
temperatures deduced from variability at 1.4 GHz are in the same range as those from 5 GHz (Table 3); this favours the intrinsic 
variability explanation, since RISS at 1.4 GHz is expected to be more pronounced than at 5 GHz due to its strong wavelength 
dependence. We therefore suggest that, although a non-negligible role played by RISS may not be ruled out in individual objects, 
especially those with relatively low amplitudes, our results for the sample as a whole should not be affected significantly. 

E. ON THE X-RAY RADIATION MECHANISMS OF RADIO-LOUD NLS1 AGN 

As one of the defining characteristics, the X-ray spectra of FSRQs are flat and extend from the soft to the hard X-ray band 
(T <~ 1.6 with a dispersion a = 0.1-0.2; see e.g. Reeves & Turner 2000), which is interpreted as inverse Compton radiation from 
the relativistic jets and beamed. In most RLQs the contribution from a hot disk corona to the observed X-rays is negligible (in 
the hard X-ray band > 2keV), except the steep-spectrum soft X-ray excess below 1 keV (e.g. Brinkmann et al. 1997; Yuan et 
al. 2000). From the data analysis above, we find a flat X-ray spectrum in J1633+4718 with Y = 1.37 ±0.49, and possibly in 
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J0849+5108 with T = 1.77^ (90% errors). Their concave optical to X-ray spectral shapes (see Figure 10) also agree with the 
inverse Compton scenario. For J0849+5 108, the LBL-like SED and flat X-ray spectrum are well expected, based on its OVV/HPQ 
properties (see AppendixB). The prominent soft X-ray excess in J1633+4718 is also typical of normal NLS1 galaxies. It is 
puzzling, however, that J1633+4718 shows a broad band SED similar to that of HFSRQs, for which a steep X-ray spectrum and 
a convex a ox -X-ray shape are generally expected. Future multi-waveband observations, especially in the X-ray band, are needed 
to confirm this peculiarity. 

As discussed above, in the presence of an excess soft X-ray component, the effective spectral slopes derived from the hardness 
ratios are generally steep. The same is expected to be true for radio-loud NLS1 galaxies, as found in J1633+4718 above. For 
reference, FSRQs at redshifts (z < 1) are found to have V ~ 2.1 with a dispersion a ~ 0.2 in the ROSAT PSPC band (Brinkmann 
et al. 1997; Yuan 1998). Therefore, we expect that a few other objects, e.g. those having FSRQ/LBL-type SEDs yet steep T 
values in the PSPC bandpass (Figure 9) would similarly have a flat underlying X-ray continuum. 

HBLs are commonly observed to show steep X-ray spectra (e.g. Worrall & Wilkes 1990; Brinkmann et al. 1996; Padovani & 
Giommi 1995; Wolter et al. 1998), which are intrinsically curved, progressively steepening with increasing photon energy (e.g. 
Perlman et al. 2005). Their X-ray emission is thought to be attributed to the high-energy tail of the synchrotron radiation, i.e. 
the synchrotron bump reaches the soft X-ray band (Ghisellini & Maraschi 1989). There are only a few flat-spectrum radio-loud 
NLS 1 AGN studied previously with multi-waveband data and, interestingly, most of these show HBL/HFSRQ-like SEDs. For 
example, RXJ 16290+4007, formerly classified as a FSRQ, was suggested by Padovani et al. (2002) to be the first confirmed 
HFSRQ for its modeled synchrotron peak at 2 x 10 16 Hz and a steep soft X-ray spectrum T ~ 3. 2MASX J0324+3410 also shows 
a continuum SED (square in Figure9; Zhou et al. 2007) and a luminosity closely resembling that of Mkn421, a well-known 
HBL; the object was claimed to be even marginally detected in TeV 7-rays (Falcone et al. 2004), as seen in some nearby bright 
HBL. In the X-ray band, 2MASX J0324+34 1 also shows a steep spectrum (T ~ 2.2, Zhou et al. 2007). 

Most of the objects falling into the HBL locus in Figure 9 show steep X-ray spectra wherever available (T = 2.0-3.3). Mean- 
while, their a ox -a x SED shapes are either convex or close to a straight line (i.e. power-law; see Figure 10). Thus, based on the 
similarity with HBL in both the SED and X-ray spectral shape, we suggest their X-ray emission to be the high-energy tail of the 
synchrotron component, as in HBL. There are a few objects lying close to the border line, known as intermediate BL Lac objects 
(IBL) where the synchrotron peaks are around 10 14 - 10 15 Hz. In them, the high-energy synchrotron tail may still reach the soft 
X-rays, as found in a few IBL objects (e.g. Tagliaferri et al. 2000; Ferrero et al. 2006). 
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TABLE 1 

Sample of very radio-loud NLS 1 galaxies 



No. 


name 


SDSS name 


Z 


u 


g 


M B 


/20cm 


/6cm 






logS 


-\/5100 a 


a\ 


fwhm(H/3) 


/(H/3 ) b 


/(OJU) b 




fflNTTV 

/V U V 












mag 


man 
mag 


mag 


ndy 


mJy 


mjy 










km s — ' 








mac 
mag 


mag 


HI 


C21 


ni 
y j j 


C41 


f51 

y j ) 


(61 


f71 
y'j 


CR1 

y°i 


C91 
y~) 


yivj 




ri2i 
y ijL ) 




T141 


M 51 


CI 61 


("171 
V 1 'J 


C1 ft} 


M91 


("2(11 


1 


J08 14+5609 


J08 1432. 11+560956.6 


0.509 


18.05 


17.89 


-24.1 


69/60 


43 


72 


-0.33 


2.53 


458 


-2.43 


2164± 59 


870±16 


54±7 


1.04 






2 


J0849+5108 


J084957.98+5 10829.0 


0.583 


19.80 


18.85 


-23.1 


344/266 


161 


202 


~0 C 


3.16 


187 


-0.20 


1811±191 


135±12 


42±3 


0.23 






3 


J0850+4626 


J085001. 17+462600.5 


0.523 


19.67 


19.02 


-22.8 


21/16 




51 




2.23 


409 


-0.86 


1251± 61 


255±10 


133± 5 


1.48 


20.80 


21.29 


4 


J0902+0443 


J090227. 16+044309.6 


0.532 


19.38 


18.92 


-23.0 


153/157 


106 




-0.30 


3.02 


520 


-0.53 


2089±205 


341 ±20 


192± 6 


1.10 


20.28 




5 


J0948+0022 


J094857.32+002225.5 


0.584 


18.59 


18.36 


-23.9 


108/70 


295 




0.82 


2.55 


510 


-1.87 


1432± 87 


341±12 


34±4 


1.22 


18.61 


19.04 


6 


J0953+2836 


J095317.09+283601.5 


0.657 


19.19 


18.94 


-23.6 


45/43 








2.71 


240 


-1.73 


2162±201 


215±14 


52±5 


1.23 






7 


J1031+4234 


J103123.73+423439.3 


0.376 


19.64 


19.36 


-21.6 


17/19 




44 




2.34 


284 


-0.55 


1642± 75 


379±14 


131± 8 


1.07 






8 


J1037+0036 


J103727.45+003635.6 


0.595 


20.04 


19.32 


-22.8 


27/28 








2.66 


225 


-1.71 


1357±107 


173± 9 


25±4 


1.58 


21.25 


22.99 


9 


J 1047+4725 


J104732.68+472532.1 


0.798 


19.16 


18.95 


-24.0 


734/789 


404 


1463 


-0.51 


3.87 


293 


-1.30 


2153±267 


421 ±32 


325± 8 


0.66 






10 


Jl 110+3653 


Jl 11005.03+365336.3 


0.630 


21.13 


20.59 


-21.8 


19/23 




45 




2.97 


74 


-1.28 


1300±339 


93±20 


32±4 


0.88 


21.35 


21.85 


11 


Jl 138+3653 


Jl 13824.54+365327.1 


0.356 


20.49 


19.52 


-21.2 


13/12 




26 




2.34 


378 




1364± 99 


393±20 


72±10 


0.69 






12 


Jl 146+3236 


Jl 14654.28+323652.3 


0.465 


19.08 


18.62 


-23.0 


15/17 




26 




2.11 


427 


-1.54 


2081 ±183 


632±39 


149±6 


0.78 






13 


J1238+3942 


J123852. 12+394227.8 


0.622 


19.89 


19.55 


-22.8 


10/10 








2.23 


175 


-1.19 


910±200 


103±18 


68±4 


1.05 






14 


J1246+0238 


J124634.65+023809.0 


0.362 


18.52 


18.31 


-22.7 


37/36 








2.38 


601 


-1.89 


1425± 68 


527±12 


115±5 


0.84 


19.01 


19.58 


15 


Jl 305+5 116 


J130522.75+51 1640.3 


0.785 


17.45 


17.28 


-25.6 


84/87 


46 


211 


-0.50 


2.34 


1264 


-1.68 


1925± 53 


1982±43 


487±13 


0.65 


17.45 


18.14 


16 


J1435+3131 


J143509.49+313147.8 


0.501 


19.87 


19.43 


-22.3 


39/43 




137 




2.87 


306 


-1.70 


1719±253 


346±36 


137±6 


1.32 






17 


J1443+4725 


J144318.56+472556.7 


0.703 


18.20 


18.11 


-24.6 


165/166 


72 


309 


-0.67 


3.07 


534 


-2.06 


1848±113 


337±18 


105±13 


1.45 


18.09 


18.57 


18 


J 1505+0326 


J150506.48+032630.8 


0.408 


18.78 


18.52 


-22.7 


365/395 


859 




0.66 


3.19 


389 


-0.83 


1082±113 


95± 8 


106±4 


1.53 






19 


J1548+3511 


J154817.92+351128.0 


0.478 


18.18 


17.91 


-23.9 


141/141 


107 


279 


-0.22 


2.84 


545 


-2.32 


2035 ± 52 


721±12 


309±7 


1.24 


18.41 


18.69 


20 


J 1633+47 18 


J163323.58+471859.0 


0.116 


17.94 


17.38 


-20.9 


63/69 


35 


107 


-0.30 


2.22 


1868 


-1.70 


909± 43 


902±21 


919± 8 


1.02 


18.20 


19.48 


21 


J 1634+4809 


J163401.94+480940.2 


0.494 


20.18 


19.53 


-22.2 


8/14 




16 




2.31 


257 


-1.60 


1609± 79 


287±11 


87± 3 


0.98 


21.03 


22.30 


22 


J1644+2619 


J164442.53+261913.2 


0.144 


18.03 


18.02 


-20.8 


88/128 


99 




-0.07 


2.65 


1788 


-1.00 


1507± 42 


1108±22 


120±4 


0.75 


18.50 


18.65 


23 


J1722+5654 


J172206.03+565451.6 


0.425 


18.66 


18.31 


-23.1 


37/43 




119 




2.37 


544 


-1.41 


1385± 36 


542±11 


323±4 


0.58 







NOTE. — Col. (2): abbreviated name of object used in this paper; Col. (3): SDSS name; Col. (4): redshift; Col. (5): SDSS u band PSF magnitude corrected for Galactic extinction; Col. (6): as col. (5), but for 
the g band; Col.(7): derived B band absolute magnitude; Col.(8): FIRST and NVSS radio flux densities at 20cm (FIRST/NVSS); Col.(9): flux density at 6 cm; Col.(10): flux density at 92cm given by WNESS; 
Col.(ll): radio spectral slope between 6cm and 20cm; Col.(12): logrifhm of radio-loudness defined as = /i,(1.4GHz)//„(4400); Col.(13): flux \f at 5100A in the rest frame; Col.(14): optical power-law slope 
(S(A) oc A"a ) ; Col.(15): line width of H/3 in FWHM; Col.(16): line flux of the broad component of H/3; Col.(17): line flux of [O III]; Col.(18): FeII-to-H/3 flux ratio (see text); Col.(19): near-UV magnitude given 
by GALEX corrected for Galactic extinction; Col.(20): the same as col.(19) but for the far-UV band. 

"Rest frame flux \f in units of 10~ 15 erg s~' cirr 2 
b Rest frame line flux in units of 10~ 17 ergs - ' cm -2 

c Given the strong variation of the 1.4GHz flux, we adopt a representing radio slope a r ~ as suggested by the radio data from Arp et al. (1979). 
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TABLE 2 

X-RAY DATA AND ANALYSIS RESULTS 



name 


X-ray source 


obs. 


r CX po(s) 

sec 


A<H.Gal 
10 20 cnT 2 


C.Rate 

10~ 2 c/s 


flux a 


r 


eV 


remark 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J0814+5609 


J081432.9+561008 


RS 


345 


4.49 


5.40 


1.71 


2.57^« 
3.36 ±6. 12 




BSC 


— 


J08 1432. 11+560956.6 


RP 


7343 


4.49 


2.04 


0.39 




WGA 


J0849+5108 


J084957.98+5 10829.0 


RP 


4472 


3.00 


1.19 


0.27 


i 77 +0.32 
-0.39 




fit 


J0850+4626 


no detection 




469 


2.64 


<4.51 


<0.98 






J0902+0443 


no detection 




362 


3.11 


<5.72 


<1.45 








J0948+0022 


J094856.9+002235 


RS 


447 


5.22 


4.10 


1.16 


2.26 ±0.64 




FSC 


J0953+2836 


no detection 




449 


1.28 


<5.52 


<0.66 








J 103 1+4234 


no detection 




498 


0.99 


<6.00 


<0.60 








J1037+0036 


J103727.45+003635.6 


RP 


10475 


4.87 


0.84 


0.16 


o 43+0.43 

49 

3.23 ±6.29 









J103727.45+003635.6 


RP 


6158 


4.87 


2.56 


0.51 






J1047+4725 


J 10473 1.4+472528 


RS 


509 


1.30 


2.13 


0.26 






FSC 


Jl 110+3653 


no detection 




315 


1.69 


<7.22 


<1.06 








Jl 138+3653 


no detection 




327 


1.75 


<6.16 


<0.94 








Jl 146+3236 


Jl 14654.0+323653 


RS 


452 


1.42 


1.67 


0.21 


9 37+0.81 
OK 




FSC 


J1238+3942 


J123852.1+394253 


RS 


508 


1.46 


5.08 


0.65 


j 4<+0.32 
^ -0.28 




BSC 


J1246+0238 


no detection 




160 


2.01 


<11.6 


<1.98 






J1305+5116 


no detection 




499 


0.98 


<5.82 


<0.58 








J1435+3131 


J143509.2+313203 


RS 


382 


1.14 


10.50 


1.15 






BSC 


J 1443+4725 


J144319.8+472541 


RS 


587 


1.47 


1.92 


0.25 






FSC 


J1505+0326 


no detection 




417 


3.89 


<6.90 


<2.14 








J1548+3511 


J154817.6+351130 


RS 


241 


2.28 


5.49 


1.05 






FSC 


J1633+4718 


J163323.58+471859.0 


RP 


3748 


1.79 


26.0 


5.44 


1.37 ±0.30 


32.5 + ^ 
30 +r2 


fit 


J1633+4718 


J163323.3+471848 


RS 


909 


1.79 


20.0 


4.72 


1 47+0.77 


fit 


J 1634+4809 


J 163400. 1+4809 15 


RS 


943 


1.64 


1.40 


0.20 






FSC 


J 1644+26 19 


J164443.2+261909 


RS 


518 


5.12 


10.6 


2.60 


2 03+ 25 




BSC 




as target 


Ch 


2946 


5.12 


18.0 b 


1.05 b 


2.19±0.17 




fit 


J1722+5654 


J172205.0+565450 


RS 


1415 


2.10 


8.11 


1.33 


2.40±0.15 




BSC 



NOTE. — Col.(l): object name; Col.(2): name of X-ray source; Col. (3): RS - the RASS, RP - ROSAT pointed observation, 
Ch - Chandra observation; Col. (4): effective exposure time in units of seconds (vegneting corrected); Col. (5): Galactic absorp- 
tion column density; Col. (6): X-ray count rate as detected with the ROSAT PSPC (except one of J1644+2619 with the Chandra 
ACIS-S; Col.(7): X-ray flux density corrected for Galactic absorption (0. 1-2.4 keV for ROSAT and 0.3-5 keV band for Chandra 
data); Col. (8): power-law photon index; Col. (9): temperature of an additional soft X-ray blackbody component; Col. (10): re- 
marks on how the X-ray data analysis is performed; fit - by spectral fitting; for the rest, the fluxes and T are derived using count 
rates and hardness ratio (see text for details) from RASS bright source catalogue (BSC), RASS faint source catalogue (FSC), 
and ROSAT source catalogues from pointed observations (RXP) and its WGA version (WGA), respectively. 

"in units of 10~ 12 ergs -1 cnT 2 

b in the 0.3-5 keV band for the Chandra ACIS-S detector 
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TABLE 3 

Radio flux variations and inferred brightness temperatures 



name b band AS AS/(S) A; log 7b ^min reference 





deg. 


GHz 


mjy 


% 


days 


K 






J0814+5609 


33.7 


1.4 


9 


14 


1273 


12.0 




2 


J0849+5108 


39.1 


1.4 


78 


25 


1257 


13.0 


2.2 


2,6 






1.4 


230 


75 


2191 


13.0 


2.2 


7,8 


J0850+4626 


39.3 


1.4 


5 


27 


1223 


11.8 




2,6 


J0902+0443 


31.2 


4.85 


10 


9 


1125 


11.1 




1,3 


J0948+0022 


38.7 


1.4 


38 


43 


1282 


12.7 


1.7 


2,6 


J1505+0326 


50.3 


4.85 


527 


72 


1125 


12.6 


1.5 


1,3 






8.3 


146 


43 


64 


14.0 


4.7 


9 


J 1548+35 11 


51.6 


4.85 


38 


44 


207 


13.0 


2.2 


1,4 


J1633+4718 


42.6 


4.85 


24 


55 


1370 


9.9 




1,5 



NOTE. — band - observing frequency; b - Galactic latitude; AS- abso- 
lute flux change between two epoches; AS/ (5)- fractional amplitudes where 
(S) is the average of the two fluxes; A?-time span of the two epoches (for 
an observation whose epoch is available accurately only to month, the maxi- 
mum length of dates is used for conservative estimation); Th-brightness tem- 
perature (logrithm); <5 m j n -minimum Doppler factor of beaming assuming an 
intrinsic maximum Tb as the inverse Compton limit 10 12 k. Reference codes: 
(1) Becker et al. (1991) (2) Becker et al. (1995) (3) Griffith et al. (1995) (4) 
Langston et al. (1990) (5) Neumann et al. (1994) (6) Condon et al. (1998a) 
(7) Willis et al. (1976) (8) Arp et al. (1979) (9) Goddard Geodetic VLBI 
GROUP'S auxiliary web pages http://lacerta.gsfc.nasa.gov/vlbi/images/ 



TABLE 4 

Broad-band effective spectral indices 3 and 
black hole masses 1 " 



name 








log(M B H/M Q ) 


J0814+5609 


0.44 


0.87 


1.75 


8.0 


J0849+5108 


0.65 


0.84 


1.23 


7.4 


J0850+4626 


0.44 


>0.68 


>1.17 


7.2 


J0902+0443 


0.60 


>0.78 


>1.16 


7.7 


J0948+0022 


0.57 


0.78 


1.22 


7.5 


J0953+2836 


0.48 


>0.74 


>1.29 


7.8 


J 103 1+4234 


0.44 


>0.71 


>1.27 


7.3 


J 1037+0036 


0.50 


0.79 


1.39 


7.3 


J 1047+4725 


0.70 


0.94 


1.45 


8.1 


Jl 110+3653 


0.56 


>0.68 


>0.92 


7.1 


Jl 138+3653 


0.47 


>0.66 


>1.07 


7.1 


Jl 146+3236 


0.39 


0.74 


1.47 


7.8 


J1238+3942 


0.41 


0.65 


1.16 


6.8 


J1246+0238 


0.41 


>0.68 


>1.25 


7.3 


J1305+5116 


0.40 


>0.78 


>1.56 


8.5 


J1435+3131 


0.52 


0.71 


1.11 


7.5 


J1443+4725 


0.51 


0.86 


1.58 


7.8 


J 1505+0326 


0.70 


>0.86 


>1.19 


6.6 


J1548+3511 


0.52 


0.80 


1.38 


7.9 


J 1633+47 18 


0.43 


0.68 


1.19 


6.3 


J 1634+4809 


0.44 


0.74 


1.35 


7.4 


J 1644+26 19 


0.51 


0.74 


1.22 


6.9 


J 1722+5654 


0.43 


0.69 


1.24 


7.4 



"Calculated using the rest frame luminosities at 
5 GHz, 2500 A, and 2 keV. The lower limits are com- 
puted for X-ray non-detections from the X-ray flux 
limits set by theRASS. 

b Virial black hole masses estimated from the width 
and luminosity of the broad H/3 lines (see text for 
details). 



